Coupled-cluster calculations for valence systems around ie O 
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We study the ground and low-lying excited states of ls O, 17 0, 15 N, and 17 F using modern two- 
body nucleon-nucleon interactions and the suitably designed variants of the ab initio equation-of- 
motion coupled-cluster theory aimed at an accurate description of systems with valence particles 
and holes. A number of properties of ls O, 17 0, 15 N, and 17 F, including ways the energies of ground 
and excited states of valence systems around 16 O change as functions of the number of nucleons, 
are correctly reproduced by the equation-of-motion coupled-cluster calculations. Within a harmonic 
oscillator basis and large effective model spaces, our results are converged for the chosen two-body 
Hamiltonians. Thus, all disagreements with experiment are, most likely, due to the degrees of 
freedom such as three-body interactions not accounted for in our effective two-body Hamiltonians. 
In particular, the calculated binding energies of ls O/ 15 N and 17 0/ 17 F enable us to rationalize the 
discrepancy between the experimental and recently published [Phys. Rev. Lett. 94, 212501 (2005)] 
equation-of-motion coupled-cluster excitation energies for the ,r = 3~ state of le O. The results 
demonstrate the feasibility of the equation-of-motion coupled-cluster methods to deal with valence 
systems around closed-shell nuclei and to provide precise results for systems beyond A — 16. 

PACS numbers: 



I. INTRODUCTION 

The way shell closures and single-particle energies 
evolve as functions of the number of nucleons is presently 
one of the greatest challenges to our understanding of 
the basic features of nuclei. The properties of single- 
particle energies and states with a strong quasi-particle 
content along an isotopic chain are moreover expected 
to be strongly influenced by the nuclear spin-orbit force. 
The latter can be retraced to contributions from both 
two-body and three-body models of the nuclear forces 
(see, for example, Refs. 0,0)- A fully microscopic ab 
initio description of masses, shell closures, excited states, 
and single-particle energies in terms of the underlying nu- 
clear forces is an unresolved problem in nuclear physics 
that awaits a satisfactory and computationally tractable 
solution. 

For light nuclei with mass numbers A < 12, both 
Green's function Monte Carlo methods and large- 
scale no-core shell-model calculations 0| provide almost 
converged benchmarks for selected two- and three-body 
Hamiltonians, where typically the models for the two- 
body nucleon-nucleon interactions reproduce the avail- 
able scattering data, while the three-body interaction 
models are normally fitted to reproduce the binding en- 
ergies of selected nuclei. The agreement with experi- 
mental data for many light nuclei is in these calcula- 
tions quite reasonable. Unfortunately, for medium-mass 
and heavier nuclei the dimensionality of the correspond- 
ing many-particle problem becomes intractable by the 
Green's function Monte Carlo methods and ab initio no- 
core shell-model techniques, and one typically has to 



resort to a simplified shell-model description within a 
smaller space, the so-called model space. In order to solve 
the corresponding many-body Schrodinger equation, one 
needs then to derive effective two and/or three-body in- 
teractions for the chosen small model space. Many-body 
perturbation theory is normally employed to derive ef- 
fective interactions but unless these interactions are 
fitted to reproduce selected properties of nuclei [(J Q, 
one cannot correctly recover the experimentally derived 
single-particle and excitation energies and shell closures 
(see, for example, Ref. @). 

Two key points make it imperative to investigate new 
theoretical methods that will allow for an accurate de- 
scription of closed- as well as open-shell nuclei with 
A ^> 12. First, present and proposed nuclear struc- 
ture research facilities will open significant territory into 
regions of medium-mass and heavier nuclei, where the 
majority of the studied nuclei will be open-shell systems 
and where many of the nuclei produced in experiment 
will be unstable or short-lived. Second, existing shell- 
model and Green's function Monte Carlo techniques have 
prohibitive computer costs that scale factorially or ex- 
ponentially with the system size. In addition to an in- 
creased dimensionality, one needs to account for the fact 
that many of the medium-mass and heavier nuclei can 
be weakly bound and couple to resonant states. More- 
over, in order to examine new nuclei that have not been 
discovered or studied before, one may not be able to rely 
on fitting the effective Hamiltonians to the experimen- 
tal data for the known nuclei, as has been tradition- 
ally done for many years. Microscopic ab initio meth- 
ods, in which nuclear properties are obtained from the 
underlying nucleon-nucleon interactions, will become in- 
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creasingly important as the new information about the 
medium-mass and heavier nuclei is obtained in various 
experiments. In addition to these practical aspects, ab 
initio calculations of nuclear properties, including, for 
example, the way the binding and excitation energies 
change as a function of the number of nucleons around 
closed-shell nuclei, may provide important new insights 
into our understanding of nuclear forces. 

Clearly, if we wish to extend ab initio methods to nu- 
clei with A ^> 12, we have to consider alternatives to the 
existing Green's function Monte Carlo and no-core shell- 
model techniques. In this work, we focus on coupled- 
cluster theory [1 OH [3 d , which is a promising can- 
didate for the development of practical methods for fully 
microscopic ab initio studies of nuclei in the A 3> 12 mass 
region. As has been demonstrated over and over again 
in numerous quantum chemis try app lications (see, e.g., 
Refs. [H Q M 13 13 13 E M,Wi 13 E3 for selected 
reviews) , coupled-cluster methods are capable of provid- 
ing a precise description of many-particle correlation ef- 
fects at the relatively low computer costs, when com- 
pared to shell-model or configuration interaction tech- 
niques aimed at similar accuracies. Based on the remark- 
able success of coupled-cluster methods in chemistry and 
molecular physics, where one has to obtain a highly accu- 
rate description of many-electron correlation effects, we 
believe that the field of nuclear physics may significantly 
benefit from the vast experience in the development of ac- 
curate and computationally efficient coupled-cluster ap- 
proximations and algorithms by quantum chemists. 

Although historicallycoupled-cluster theory originated 
in nuclear physics [fj, I3l > its applications to the nuclear 
many-bod y p roblem have been relatively rare (see, e.g., 
Refs. [24|, and references therein), particularly 

when compared to quantum chemistry. For many years, 
part of the problem has been an inadequate understand- 
ing of nucleon-nucleon interactions and lack of adequate 
computer resources in the 1970s and 1980s. This situ- 
ation has changed only in the last few years. The suc- 
cessful construction of realistic nucleon-nucleon poten- 
tials (cf., e.g., Refs. |3 13 13 HI) and spectacular 
improvements in computer technology have led to re- 
newed interest in applying coupled-cluster methods in 
ab initio nuclear physics calculations. In particular, us- 
ing bare interactions, Mihaila and Heisenberg performed 
impressive coupled-cluster calculations for the binding 
energy and the electron scattering form factor of 16 O 
[3ll |3 13 13] ■ We have taken an alternative route and 
combined a few basic coupled-cluster techniques, devel- 
oped earlier in the context of electronic structure studies 
by quantum chemists, with the renormalized form of the 
Hamiltonian to determin e g round and selected excited 
states of 4 He and 16 O |^|^|^|^|3!4|4(j, demonstrat- 
ing promising results when compared with the results of 
the exact shell-model diagonalization in the same model 
space [3^ | and, at least for some properties, with the ex- 
perimental data H^, H3, 13, 13 113 • In particular, in our 
most recent study of the ground and excited states of 



16 O [33, we obtained fully converged results which are 
very close to those obtained with more expensive large- 
scale no-core shell- model calculations of Navratil and col- 
laborators ^lj. This has been possible thanks to the 
use of the elegant diagram factorization techniques de- 
veloped by quantum chemists 01 > which lead to almost 
perfectly vectorized [3EUEI an d highly scalable paral- 
lel 35] computer codes, enabling routine calculations for 
systems in the A ~ 20 region with large single-particle 
basis sets, including seven or even eight major oscillator 
shells (336 and 480 single-particle states, respectively). 
It should be emphasized that although we are still in the 
early stages of developing a library of efficient, general- 
purpose coupled-cluster programs for nuclear structure 
applications, there are already several differences be- 
tween our approach to nuclear coupled-cluster calcula- 
tions and th e appro ach pursued by Mihaila and Heisen- 
berg |3 13 03 13| First of all, Mihaila and Heisenberg 
used bare interactions, making the convergence with the 
number of single-particle basis states very slow, whereas 
we use the renormalized form of the Hamiltonian exploit- 
ing, for example, a no-core G-matrix theory |35j . which 
leads to a rapid convergence of binding and excitation 
energies and other nuclear properties with the number of 
major oscillator shells in a basis set |35l l37l l38l I39tl40 | . 
Second, we are able to calculate ground as well as excited 
states, not just the ground-state properties of closed-shell 
nuclei examined by Mihaila and Heisenberg. Finally, 
as mentioned above, our coupled-cluster computer codes 
have been developed using diagram factorization tech- 
niques, which minimize the CPU operation count, rather 
than the commutator expansions used by Mihaila and 
Heisenberg and the Bochum school. The coding style 
adopted by us is similar to that practiced by the lead 
developers of coupled-cluster methods in chemistry. In 
particular, we put an emphasis on the general-purpose 
character of our codes, meaning that the only essential 
input variables are the number of particles and the ma- 
trix elements of the Hamiltonian in some single-particle 
basis set. 

Our initial coupled-cluster calculations [3 13 13 13 
have focused on closed-shell nuclei. However, the 
long-term objective of our research program is to study 
open-shell nuclei with one or more valence nucleons. We 
would like, for example, to examine how the binding and 
excitation energies vary with the number of nucleons in 
valence systems around closed-shell nuclei. This is par- 
ticularly interesting, when we examine the A = 15 and 
A = 17 nuclei around 16 0. For example, the splittings 
between the (3/2)^ and (1/2)^ states in 15 O and 15 N 
and the splittings between the (3/2)^ and (5/2)^ states 
in 17 and 17 F should arise from the nuclear spin-orbit 
force, which may or may not be affected by three-nucleon 
interactions. One would like to examine such issues by 
comparing the results of converged ab initio calculations 
employing two-body interactions with the experimental 
energy spacings. This requires, however, an appropri- 
ate extension of the usual single-reference ground-state 
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coupled-cluster theory d 0, ITTl IT^ to ground and ex- 
cited states of valence systems around closed-shell nuclei. 

In this paper, we examine, for the first time, the ap- 
plicability of two quantum-chemistry-inspired coupled- 
cluster approaches, referred to as the particle-attached 
(P A) ( in chemistry, electron-attached or EA [45t l46l 
I47L l49| l and p article-remove d ( PR) (in chemistry, 
ionized or IP [lH |H |H El El El El El Eif ) 
equation-of- motion coupled-cluster (EOMCC) methods 
[45L El E3i m the converged calculations of the bind- 
ing and excitation energies of the j4 = 15 ( 15 0, 15 N) 
and A — 17 ( 17 0, 17 F) nuclei. For these calculations, 
we use modern nucleon-nucleon interactions derived from 
the effective-field theory jH E3, such as N 3 LO 
and their slightly older phenomenological counterparts, 
including the charge-dependent Bonn interaction model 
(CD-Bonn) |2S[ and the Vis model of the Argonne group 
|27|. In the PA- and PR-EOMCC methods, one cal- 
culates ground and excited states of the (A + 1)- and 
(A — l)-particle systems by diagonalizing the similarity- 
transformed Hamiltonian of the coupled-cluster theory, 
resulting from the ground-state calculations for the A- 
particle closed-shell nucleus in the relevant (A + 1)- and 
(A— l)-particle subspaces of the Fock space. As shown in 
this paper, the PA- and PR-EOMCC approaches provide 
us with practical computational techniques for poten- 
tially accurate ab initio studies of valence systems around 
the closed-shell nuclei that may provide several impor- 
tant insights into the effects of the underlying nucleon- 
nucleon interactions on the calculated properties of such 
systems. In addition to the converged PA- and PR- 
EOMCC results for the A = 15 and A = 17 nuclei ob- 
tained with three different types of nucleon-nucleon inter- 
actions (N 3 LO, CD-Bonn, and Vis), we provide several 
details of the PA-EOMCC, PR-EOMCC, and underlying 
ground-state coupled-cluster calculations, including the 
factorized forms of the relevant amplitude equations that 
lead to highly efficient, fully vectorized computer codes 
applicable to large single-particle basis sets and masses 
in the A ~ 20 - 40 region. 

This paper is divided into four sections. In Sec. [H] 
we present our formalism for deriving an effective two- 
body Hamiltonian for coupled-cluster calculations, which 
takes into account short-range nucleon-nucleon correla- 
tions, and present the details of the PA-EOMCC and 
PR-EOMCC theories that enable us to deal with va- 
lence systems around closed-shell nuclei within the frame- 
work of the single-reference coupled-cluster formalism. 
The results of PA-EOMCC and PR-EOMCC calculations 
for the 15 0, 17 0, 15 N, and 17 F nuclei are discussed in 
Sec. 11111 and the conclusions and perspectives are out- 
lined in Sec. US The factorized forms of the PA-EOMCC 
and PR-EOMCC equations for the (^4 ± l)-particle sys- 
tems and the corresponding ground-state coupled-cluster 
equations, exploited in this work, are given in the Ap- 
pendix. 



II. COUPLED-CLUSTER EQUATIONS FOR 
VALENCE SYSTEMS 

This section serves the aim of presenting the coupled- 
cluster theories for open-shell nuclei used in this work, 
with an emphasis on the PA-EOMCC and PR-EOMCC 
methods mentioned in the Introduction. These meth- 
ods are designed to handle systems with one valence 
particle or one valence hole. Since our effective model 
spaces for the coupled-cluster calculations, which enable 
us to obtain converged PA-EOMCC and PR-EOMCC re- 
sults, involve up to eight major harmonic oscillator shells 
(480 uncoupled single-particle basis states) and 15 to 17 
explicitly correlated nucleons, we focus on the compu- 
tationally efficient formulation of the PA-EOMCC and 
PR-EOMCC methods that makes such large-scale nu- 
clear structure calculations manageable. 

Our theoretical considerations start with the introduc- 
tion of an appropriate two-body effective interaction for 
the large-scale coupled-cluster calculations. This is be- 
cause the majority of modern nucleon-nucleon interac- 
tions, including the N 3 LO, CD-Bonn, and Vis potentials 
examined in this study, include repulsive cores that would 
require calculations in extremel y la rge model spaces to 
reach converged results 0, El El E3 ■ ^ n or der to re- 
move the hard-core part of the interaction from the prob- 
lem and allow for realistic calculations in manageable 
model spaces with seven or eight major oscillator shells, 
one has to renormalize the bare interactions. The rele- 
vant information about the method used by us to renor- 
malize the N 3 LO, CD-Bonn, and Vis Hamiltonians and 
to generate the final effective Hamiltonians corrected for 
the center-of-mass contaminations, which can be used in 
the PA-EOMCC, PR-EOMCC, an d oth er coupled-cluster 
calculations, are discussed in Sec. Ill Al 

As mentioned in the Introduction, the PA- and PR- 
EOMCC methods for valence systems are based on an 
idea of diagonalizing the similarity-transformed Hamil- 
tonian of the coupled-cluster theory, resulting from the 
ground-state calculations for the A-particle closed-shell 
nucleus, in the relevant (A+ 1)- and (A~ l)-particle sub- 
spaces of the Fock space. Thus, before introducing the 
details of the PA- and PR-EOMCC calculations for the 
(A + 1)- and (A — l)-particle systems, we provide the 
most essential information about the underlying closed- 
shell coupled-cluster calculations that precede the PA- 
and PR-EOMCC steps. This is done in Scc. lHBl 

Finally, in Sec. Ill CI we discuss the most essential de- 
tails of the PA- and PR-EOMCC calculations for the 
ground- and excited states of the (A + 1)- and (A — 
l)-particle valence systems around the closed-shell A- 
particle nucleus. The final working equations for the clus- 
ter and excitation amplitudes, which define the coupled- 
cluster, PA-EOMCC, and PR-EOMCC approximations 
implemented in this work and which lead to highly effi- 
cient computer codes, are shown in the Appendix. 
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A. Effective Two-Body Interaction for 
Coupled-Cluster Calculations 

In the PR-EOMCC calculations for 15 and 15 N, the 
PA-EOMCC calculations for 17 O and 17 F, and the un- 
derlying closed-shell coupled-cluster calculations for the 
ground-state of 16 which precede the PR-EOMCC and 
PA-EOMCC calculations, we used the following three 
nucleon-nucleon interactions: N 3 LO jsjj, CD-Bonn [28| . 
and Via |27j . The Coulomb interaction was included in 
all of our calculations (to distinguish between 15 0/ 17 
and 15 N/ 17 F). In order to remove the hard-core part of 
the interaction, which would require calculations in ex- 
tremely large model spaces consisting of dozens of ma- 
jor oscillator shells to reach reasonably converged re- 
sults [^J |U H3, E3, an d enable realistic calculations 
in manageable model spaces, we follow the procedure 
exploited in our earlier work |H HI El HI HI E3- 
Thus, we renormalize the Hamiltonian through a no- 
core G-matrix procedure, described in considerable de- 
tail in Refs. 0, |35|. The no-core G-matrix approach 
introduces a starting-energy (u>) dependence in the effec- 
tive two-body matrix elements G(w) defining the renor- 
malized two-body interactions (obtained by analyzing 
the exactly solvable proton-proton, proton-neutron, and 
neutron-neutron two-body problems), but much of the 
u) dependence can be eliminated through the use of the 
Bethe-Brandow-Petschek theorem |(3lJ and the appropri- 
ate summation of the class of folded diagrams to infinite 
order at a given starting energy (see Refs. 0, |3j| for fur- 
ther information). For nuclei like 16 0, the dependence on 
the chosen starting energy u> is weak (almost none when 
seven or eight major oscillator shells are employed |35|). 
It introduces an uncertainty of 0.1 — 0.2 MeV per particle 
for the binding energies. 

After renormalizing bare interactions with the G- 
matrix approach, our effective Hamiltonian is given by 
the formula 

fleffH = flo + G(fa>), (1) 

where Ho is the total kinetic energy of the nucleons. To 
complete the process of preparing the Hamiltonian for 
coupled-cluster calculations, we correct the renormalized 
Hamiltonian H e s(u>), Eq. Q, resulting from exploiting 
the no-core G-matrix procedure, for center-of-mass con- 
taminations using the expression 

H = H(u), Pcom) = H c g(u}) + /3coM#CoM 

= Za a ° a P + 3 v^pcfaPasa^, (2) 

where z@ = (a\z\(3) and y^g = ( a f3\ v \jd) ~ ( a P\ v \fij) 
are the relevant one- and two-body matrix elements in a 
single-particle basis set {\a)} and a a (a a ) are the usual 
creation (annihilation) operators. Here and elsewhere in 
the present paper, we use the Einstein summation con- 
vention over repeated upper and lower indices. The pa- 
rameter PcoM is chosen such that the expectation value of 



the center-of-mass Hamiltonian HcoM with the ground- 
state coupled-cluster wave function, (Hc m), obtained 
for the /?coM-dependent Hamiltonian H, Eq. J2J, is 0.0 
MeV. This can be done by relying, for example, on the 
Hellmann-Feynman theorem and calculating {Hc m) as 
the first derivative of the coupled-cluster energy with re- 
spect to PcoM- Once we know the values of (-Hcom) at 
various /?coM values, we can easily identify the optimum 
/3coM value at which {Hcom) becomes 0.0 MeV [6(J. As 
pointed out in our earlier papers jH H^, ^3 , one of the 
advantages of this procedure is the ease of separation 
of intrinsic and center-of-mass contaminated states by 
analyzing the dependence of coupled-cluster energies on 
/3coM- As shown in Refs. (23>E3|> the physical states ob- 
tained in coupled-cluster calculations are virtually inde- 
pendent of /3coM 7 while the center-of-mass contaminated 
states show a strong, nearly linear dependence of excita- 
tion energies on flcoM- 

We are currently working on alternative approaches 
to the effective interaction. In one of these alternatives, 
instead of relying on the G-matrix method, we will con- 
struct the renormalized Hamiltonian for coupled-cluster 
calculations with the help of the Lee-Suzuki procedure 
[62l |63L|64 | . exploited in the no-core shell- model ap- 
proach |J|. This procedure will eliminate the starting- 
energy dependence from our calculations. In particular, 
we will investigate the differences between our G-matrix 
approach and the no-core approach based on the Lee- 
Suzuki transformation in the forthcoming work. Both 
methods have the appealing feature that the effective in- 
teraction is properly renormalized as the size of the har- 
monic oscillator basis is increased, approaching the bare 
Hamiltonian in the limit of an infinite single-particle basis 
set. We will also study coupled-cluster applications us- 
ing the the Viowk effective interaction approach (65|. In 
this method, one uses the cutoff in momentum space to 
soften the nucleon-nucleon interaction. The Lee-Suzuki- 
based and G-matrix-based no-core procedures rely on the 
harmonic oscillator basis cutoff only, while the Viowk pro- 
cedure uses the momentum-space cutoff in addition to the 
basis set cutoff. The three methods will produce identi- 
cal results in the limit of an infinite basis set and infinite 
momenta. 

Once the one- and two-body matrix elements of the 
center-of-mass-corrected renormalized Hamiltonian H, 
Eq. (J2J, are determined, we solve the nuclear many- 
body problem using coupled-cluster theory. In order to 
construct the coupled-cluster equations for the closed- 
shell A-body system and the related PA-EOMCC and 
PR-EOMCC equations for the (A+ 1)- and (A - l)-body 
nuclei in the computationally most efficient way, as dic- 
tated by the factorized form of these equations discussed 
in Secs. lllBllII CI and the Appendix, we sort the one- and 
two-body matrix elements of H according to the particle- 
hole (p-h) character of the single-particle indices that la- 
bel them prior to the coupled-cluster work. For example, 
the two-body matrix elements u^t defining H, Eq. J3J), 
are sorted out into six groups corresponding to the fol- 
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lowing p-h types of single-particle indices a, /3, 7, <5: hhhh, 
hhhp, hhpp, hphp, hppp, pppp. The same philosophy of 
sorting matrix elements according to the p-h character of 
the relevant single-particle indices is applied to one- and 
two-body matrix elements of the similarity-transformed 
Hamiltonian of coupled-cluster theory discussed in the 
next subsection. This is a common practice in the most 
efficient implementations of coupled-cluster methods by 
quantum chemists and we follow the same recipe here. 

B. Brief Synopsis of the Single- Reference 
Coupled-Cluster Theory and the Basic CCSD 
Approximation 

As mentioned earlier, the PA-EOMCC and PR- 
EOMCC calculations for the (A + 1)- and (A - 1)- 
particle valence systems rely on the diagonalization of 
the similarity-transformed Hamiltonian, obtained in the 
single-reference coupled-cluster calculations for the A- 
particle closed-shell nucleus, in appropriate subspaces of 
the Fock space. Thus, before introducing the PA- and 
PR-EOMCC methods for the (A+l)- and (A-l)-particle 
systems, we provide the most essential information about 
the underlying closed-shell coupled-cluster calculations 
that precede the PA- and PR-EOMCC steps. 

The single-reference coupled-cluster theory 0, |n], 
is based on the exponential ansatz for the ground- 
state wave function of the ^4-body system, 

l^> = e r( », (3) 

where T^ A > is the cluster operator (a p-h excitation op- 
erator) and |$) is the corresponding reference determi- 
nant (defining the Fermi vacuum) obtained by perform- 
ing some mean-field calculation or by simply filling A 
lowest-energy single-particle states (this is what we have 
done in the calculations discussed in this paper). Here 
and elsewhere in the present paper, we use superscripts, 
such as (^4), which indicate the number of particles in 
a system under consideration, at the relevant operators 
and energies. Normally, when A remains fixed through- 
out the entire calculation, this is not essential, but in this 
paper we deal with systems with different mass numbers 
(the ^4-body as well as the (A + l)- and (A — l)-body 
systems), so that it is useful to indicate the number of 
particles in a many-body system of interest at the most 
essential mathematical quantities to avoid confusion. 

Formally, Eq. © is a direct consequence of the 
connected-cluster theorem, first clearly stated by Hub- 
bard |66j . which is, in turn, related to the linked clus- 
ter theorem of many-body perturbation theory |66l l6?l 
l68l l6^| . According to the connected-cluster theorem, 
the cluster operator T( A ) generates connected wave func- 
tion diagrams summed through infinite order. Opera- 
tionally, T {A *> is a simple many-body excitation operator, 
which in all standard coupled-cluster approximations is 
truncated at a given (usually low) p-h excitation level 
M < A. An example of the standard coupled-cluster 



method is the CCSD (coupled-cluster singles and dou- 
bles) approach [Tol FfH W\ , which is used in this work 
to obtain the ground-state information for the closed- 
shell 16 O nucleus We label coupled-cluster methods by 
the standard acronyms adopted by chemists who have 
led the development of coupled-cluster approaches for 
over 30 years now and they are, by far, the most fre- 
quent users of coupled-cluster approaches; see, for exam- 
ple, Refs. [T3.l25| for the relevant historical remarks. In 
this case, M = 2 and the cluster operator is approx- 
imated by 

T {A) (CCSD) = T (A) (2) = Ti + T 2 , (4) 

where 

T x = f a a a ai (5) 

and 

T2 = \t%a a a b a jai (6) 

are lp-lh or singly excited and 2p-2h or doubly excited 
cluster components, t l a and t^ b are the corresponding 
singly and doubly excited cluster amplitudes, and 
(a, b, . . .) are the single-particle states occupied (unoccu- 
pied) in the reference determinant |<£>). The general form 
of the truncated cluster operator, defining a standard 
single-reference coupled-cluster approximation character- 
ized by the excitation level M, is 

M 

T^(M) = Y / T n , (7) 

n=l 

where 

Tn= 00 "< •••"<: (8) 

(n — 1, . . . , M ) are the many-body components of 
T^ A \M), and are the corresponding cluster am- 

plitudes. 

The cluster amplitudes t^ 1- are determined by 
solving a coupled system of nonlinear and energy- 
independent algebraic equations of the form: 

($%-£\H N (M)\$) = 0, *!<■■•< i n , a x < ■ ■ ■ < a n , 

(9) 

where n = 1, . . . , M, 

H N (M)=e~ T<A, ^H N e TiA> W = (H N e ^ ) W) c 

(10) 

is the similarity-transformed Hamiltonian of the coupled- 
cluster theory truncated at Mp-Mh excitations, sub- 
script C designates the connected part of the cor- 
responding operator expression, and = 
a a i . . . a an ai n ■ ■ ■ ai t |<£>) are the np-nh or n-tuply excited 
determinants relative to |<&). The operator Hjq entering 
Eq. (|10fl is the Hamiltonian in the normal-ordered form 
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relative to the ^-particle Fermi vacuum reference state 

i*>, 

Hn = H — = fi N[a a a p ] + \ N[a a a^a s a y ], 

(11) 

where f@ = (a\f\{3) = zP+v^l are matrix elements of the 
Fock matrix and N[- ■ ■] designates the normal product. 
In particular, the standard CCSD equations for the singly 
and doubly excited cluster amplitudes t % a and t^ b , defining 
Xi and T 2 , respectively, can be written as 

($? (Fat (CCSD) | $) = 0, (12) 
<$#|fl*(CCSD)|$> = 0, i<j,a<b, (13) 

where 

Jj^CCSD) = i/ Ar (2) = e - T( " ) ( CCSD ) J ff Jv e T(A, ( CCSD ) 
= (H N e T ^ ccs ^) c (14) 

is the similarity-transformed Hamiltonian of the CCSD 
approach. As mentioned in Sec. Ill Al we do not use 
the bare Hamiltonian in our nuclear structure calcula- 
tions. Thus, the Hamiltonian H used to construct the 
similarity-transformed Hamiltonian Hm{GCSD) for the 
nuclear structure calculations discussed in this paper is 
replaced by the renormalized form of the Hamiltonian, 
Eq. J2J, resulting from the no-core G-matrix calcula- 
tions. 

The system of coupled-cluster equations, Eq. 10, is 
obtained in the following way (suggested by Cfzek 11]). 
We first insert the coupled-cluster wave function l*^), 
Eq. ©, into the A-body Schrodinger equation, 

H N \*W)=AEW\¥ A \ (15) 

where 

AE { Q A) = E (A) - (16) 

is the corresponding energy relative to the reference en- 
ergy {^\H\^) 1 and premultiply both sides of Eq. (|15fl on 

the left by e _T< ' to obtain the connected-cluster form 
of the Schrodinger equation 0, 0, 0, |7^| , 

H n \$)=AeW\$), (17) 

where 

H N = e- T( H N e Ti = {H N e Ti > (18) 

is the similarity-transformed Hamiltonian. Next, we 
project Eq. I|17fl . m which is replaced by its ap- 

proximate form T^(M) : Eq. J7J), onto the excited de- 
terminants |$"* "'""), with n = 1, . . . , M, corresponding 
to the p-h excitations included in T^ A \M). The excited 
determinants l^?" 1 '""") are orthogonal to the reference 
determinant |$), so that we end up with the nonlinear 



and energy-independent algebraic equations of the form 
of Eq. ©. Once the system of equations, Eq. @, is 
solved for T< A )(M) or d'.'.'X ( or Ec l s - 03 and 
are solved for T\ and T 2 or t l a and t^ b ), the ground-state 
coupled-cluster energy is calculated using the equation 

E [A) {M) = ($|2T|$) + AE (A) {M) 

= ($|#|$) + ($|Fjv(M)|$) 

= (<Z>\H\$) + ($\H N>close (M)\$), (19) 

where Hn, c \osc(M) is the closed part of Hn(M) which 

is represented by those diagrams contributing to 
that have no external (uncontracted) Fermion lines (as 
opposed to the open part of Bn(M) which is represented 
by the diagrams having external or uncontracted Fermion 
lines; cf. Sec. Ill C|) . It can easily be shown that if H (in 
our case, the renormalized Hamiltonian defined by Eq. 
Pjl) does not contain higher-than-two-body interactions 
and 2 < M < A, we can write 

E [A) {M) = ($| J ff|$) + ($|[ J ff Ar (T 1 +T 2 +iT 1 2 )] c |$). (20) 

In other words, we only need T% and T 2 clusters to cal- 
culate the ground-state energy E [A) (M) of the A-body 
(A > 2) system even if we solve for other cluster compo- 
nents T n with n > 2. Equation l)19|l can be obtained by 
projecting the connected-cluster form of the Schrodinger 
equation, Eq. 117fl . on the reference configuration |i>) 
and replacing T^ A > by T^ A \M). In fact, the nonlinear 
character of the system of coupled-cluster equations of 
the form of Eq. © does not mean that the resulting 
equations contain very high powers of T^ a \M). For ex- 
ample, if the Hamiltonian H (in our case, the renormal- 
ized Hamiltonian obtained using the G-matrix technique) 
does not contain highcr-than-pairwise interactions, the 
CCSD equations for the Ti and T 2 clusters, or for the 
amplitudes t l a and t l ^ b that represent these clusters, be- 
come 

($?\[H N (1+T 1 +T 2 + ±T?+T 1 T 2 + ±T?)} C \$) = 0, (21) 

{$$\[H N (\ +T 1 +T 2 + \Tl + T{T 2 + \Tl 

+3^ + l^Ta + £l?)] c |$) = Q. (22) 

In general, if the Hamiltonian H contains two-body inter- 
actions only, there are no terms in cluster components T n 
that are higher than quartic terms in the coupled-cluster 
system, Eq. ©, independent of the truncation scheme 
M used to define T^(M). This is a purely mathemati- 
cal statement, resulting from the fact that we must con- 
nect external lines of Hjy with the vertices representing 
the many-body components of the cluster operator 
to determine the connected operator product represented 
by if/v, Eq. (XBt ; and not the result of some arbitrary 
truncation of the exponential coupled-cluster wave func- 
tion. 
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In this work, we apply the CCSD approach to the 
closed-shell 16 nucleus and, what is even more impor- 
tant in this particular study, to obtain the correspond- 
ing similarity-transformed Hamiltonian _ffjv(CCSD), Eq. 
(|T4*|) . which is used in the subsequent PR-EOMCC and 
PA-EOMCC calculations for the ground and excited 
states of the 15- and 17-particle nuclei. As shown in our 

earlier papers [33,|3a|33,|33,|33,|40|, ^ e CCSD method is 
perfectly adequate for a highly accurate description of the 
ground state of 16 O, providing the converged description 
of the ground state of 16 O for the Hamiltonians contain- 
ing two-body interactions (the three-body clusters T3 are 
negligible [3a, 133, lMJ Efl ) • The highly accurate descrip- 
tion of the ground state of 16 O is essential for obtaining 
a well-balanced description of the 15- and 17-particle va- 
lence systems around 16 O. The ability of the CCSD ap- 
proach to provide a highly accurate description of 16 O is 
also useful from the practical point of view. When prop- 
erly implemented, the CCSD approach is characterized 
by the relatively inexpensive ri^n^ steps, where n and 
n u are the numbers of occupied and unoccupied orbitals, 
respectively, in the single-particle basis set, making it ap- 
plicable to systems with the larger values of A. Our codes 
are already efficient enough to perform the CCSD calcu- 
lations with A ~ 20 — 40 and up to eight major oscillator 
shells, which at least for A = 16 is sufficient to obtain 
the converged description 

The explicit and computationally efficient form of 
the CCSD and other coupled-cluster equations that can 
be used in routine calculations for many-body systems 
with larger A values and larger basis sets, in terms of 
one- and two-body matrix elements of the (renormal- 
ized form of the) Hamiltonian and cluster amplitudes 
*oV "*o ( m tne CCSD case, t\ and i^ b ) can be most conve- 
niently derived by applying diagrammatic techniques of 
many-body theory combined with diagram factorization 
methods which yield highly vectorized computer codes 
[23L l35l l42l I43L Ej . Once these equations are properly 
coded , w e solve them using iterative procedures, such as 
DIIS |l3 (see Refs. 0|za). The explicitly connected 
form of the coupled-cluster equations, such as Eqs. © 
or 121fl and l|22[l . guarantees that the process of solving 
these equations leads to connected terms in cluster com- 
ponents T n and connected terms in the energy Eq A ^ (M) , 
independent of the truncation scheme M used to define 
T { AS > (M). The absence of disconnected terms in (M) 
and E^ A \m) is essential to obtain the rigorously size- 
extensive results. The computationally efficient form of 
the CCSD equations for the case of pairwise interactions, 
which can be applied in large-scale nuclear structure cal- 
culations such as those discussed in this work, in terms 
of one- and two-body matrix elements of the similarity- 
transformed Hamiltonian JJtv(CCSD) that serve as the 
most natural intermediates for setting up these equations 
and other recursively generated intermediate quantities 
that are used to obtain a perfectly vectorized computer 
code, is given in the Appendix. 



C. Equation-of-Motion Coupled-Cluster Methods 
for Valence Systems: The PA-EOMCCSD and 
PR-EOMCCSD Approximations and their 
Implementation for Nuclear Structure Calculations 

In addition to providing natural intermediates for 
setting up coupled-cluster equations, the use of 
the similarity-transformed Hamiltonians, Hn(M) or 
H N (CCSD), Eqs. JTJH OT {HJ, respectively, in coupled- 
cluster calculations provides a natural mechanism for ex- 
tending the ground-state coupled-cluster theory to ex- 
cited states of a given A-body system or to ground and 
excited states of the (A + k)- or (A — fc)-particle sys- 
tems obtained by attaching k particles to or removing k 
particles from the A-particle closed-shell core. This can 
be most efficiently done by exploiting the EOMCC for- 
malism 0, 15IH [57l a nd its PA-EOMCC (in chemistry, 
EA-EOMCC glil ElEIJiir) and PR-EOMCC (in 
chemistry, IP-EOMCClllMll HI El IH |H E3 ) 
variants, and their various multiply attached and multi- 
ply removed or ionized (cf., e.g., Refs. l-ill l7rj I??], ITsf ) 
extensions (see, also, Refs. |79l I80T l8ll l82t l83j and 
HI IH HE El El El] for the related linear response 
and symmetry-adapted cluster configuration interaction 
formalisms, respectively). In all of these methods, we ob- 
tain excited states l 1 ^^) (/U > 0) of the A-particle sys- 
tem or ground and excited states \^ A±k ^} (fi > 0) of the 
A ± fc-particle (k > 0) systems by applying the suitably 
defined excitation (R A ) or particle-attaching/particle- 
removing (Rff^ ) operator to the ground state obtained 
in the single-reference coupled-cluster calculations for the 
closed-shell A-body system. Operators and i?^ A±fe ^ 
are obtained by diagonalizing the similarity-transformed 
Hamiltonians, such as Hn(M), Eq. H10|) . in the case 
coupled-cluster theory truncated at M-tuple excitations, 
or i?jv(CCSD), Eq. O, in the CCSD case, in the rel- 
evant A-particle and (A ± fc)-particle subspaces of the 
Fock space. For example, in the EOMCCSD approxima- 
tion [56ll57| . which is a basic EOMCC approximation for 
the calculations of excited states of the A-particle system, 
we represent excited states \E')^ 4 ' ) ) as 

|^) = 4^|^ A) )=<)e^|$) (23) 

and replace by the cluster operator T^(CCSD), 

Eq. Q, obtained in the CCSD calculations, and R^ by 

4^(CCSD) ee 4 A >(2) = R^ Q + R P ,i + i^,2, (24) 
where 

R P ,o = r 1, (25) 

R pA = r* a a a i: (26) 

and 

Rp, 2 = lr%a a a h a jai (27) 
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are the reference, lp-lh, and 2p-2h components of 
R (A) (GCSD), and r , r*, and r% are the corresponding 
excitation amplitudes (1 in Eq. H25[l is a unit opera- 
tor). In a more general case of coupled-cluster theory 
truncated at M-tuple excitations, where is approx- 
imated by T( A \M), Eq. Q, the corresponding p-h ex- 
citation operator Rff^ in Eq. i|23|) is approximated by 


M 

Ri A) (M) = J2^,n, (28) 
where, in analogy to T„, 

= (^) <V.V.a„ a ai ' ' ' a a ^„ • • • a h , (29) 

for n > 1, and R^q is defined by Eq. I|25|) . In the 
EOMCC methods for excited states, sometimes referred 
to as the EE-EOMCC (excitation energy EOMCC) ap- 
proaches the p-h excitation amplitudes r* 1 ';;' 1 ™ defin- 
ing R^\m) and the corresponding excitation energies 
lu (A) (M) = E (A) (M) - E {A) {M) of the A-body system 
of interest are obtained by diagonalizing the similarity- 
transformed Hamiltonian H^(M) or, more precisely, its 
open part 

H N , opcn (M) = {H N e T(A)(M ^) c ^ ca 

= e- T(A) WH N e T(A) W-H N , close (M) 
= e -T (A HM) HNe T^(M) _ AE { A) (M), 

(30) 

which has at least two external Fermion lines, in the sub- 
space of the A-particle Hilbert space spanned by the same 
excited determinants l^ 1 '""™) with n = 1,...,M that 
are used to solve the underlying ground-state coupled- 
cluster calculations (AE^ A) (M) in Eq. (J3UJI is the 
ground-state coupled-cluster energy relative to the ref- 
erence energy cf. Eq. CHI)- In particular, 
the r l a and r^ b amplitudes of the standard EOMCCSD 
(or EE-EOMCCSD) theory and the corresponding ex- 
citation energies w^^CCSD) of the A-body system are 
obtained by diagonalizing the open part of the similarity- 
transformed Hamiltonian of the CCSD approach, 

#7V,opcn(CCSD) = (i/A r e T(A)(CCSD) ) C ,opcn, (31) 

in the subspace spanned by the singly and doubly excited 
determinants \§f) and 1$"^) used to set up and solve the 
ground-state CCSD equations. We used the EOMCCSD 
approach to calculate the selected excited states of the 
4 He and 16 nuclei in Refs. |M El S |H H3 . 

The idea of diagonalizing the similarity-transformed 
Hamiltonians H Ny0pcn (M), Eq. J3DJ, or H N , opcn (CCSD) , 
Eq. can be extended to ground and excited states of 
open-shell nuclei with (A ± k) particles by replacing the 



particle-conserving p-h excitation operator R^ in Eq. 
(|23[1 by the suitably defined particle-attaching or particle- 
removing operator i?)/^ . For example, in the basic 
PA-EOMCCSD (in chemistry, EA-EOMCCSD) jH Ejl 
and PR-EOMCCSD (in chemistry, IP-EOMCCSD) [H 
HE EHH approaches (cf., also, Ref. HI, m exploited 
in this work, we define the wave functions of the (A + 1)- 
and (A — l)-particle systems, respectively, as 

H A± V) = Ri A± Ve TW \Z), (32) 

where is approximated by T( A )(CCSD), Eq. gj, ob- 
tained in the CCSD calculations for the A-particle closed- 
shell system, and i?^ 4+1 ^ and Rff -1 ' are replaced by the 
appropriately truncated operators, 

R {A+1 \2p-lh) = R^ lp + R^. 2p - lh = r a a a + \rj b a a a b a 3 

(33) 

and 

R {A - 1) (2h-lp) = R^ lh + R^h-lp = r l ai + ±r lj b a b a 3 a u 

(34) 

respectively, which generate the (A + 1)- and (A — 1)- 
particle states from the A-particle CCSD wave function 
e Tl+T2 |<I>). The lp and 2p-\h amplitudes r a and r i, 
respectively, entering Eq. (|33|l and defining the PA- 
EOMCCSD model, and the Ih and 2h-\p amplitudes 
r l and r lJ b , respectively, entering Eq. (|34|l and defining 
the PR-EOMCCSD model, are determined by solving the 
eigenvalue problem 

(H N , opcn R^) c m = J A± VR^\t>), (35) 

where -/?7v,open is replaced by the similarity-transformed 
Hamiltonian of the CCSD theory, £f/v,o P en(CCSD), Eq. 
(|31|l . in the relevant subspaces of the (A + l)- and (A—l)- 
particle subspaces, J^^ A+1 ^ and J^' 4-1 ), respectively, of 
the Fock space. The subspace of ,J4?( A+1 ) used to solve 
the PA-EOMCCSD eigenvalue problem is spanned by the 
$ a ) = a «|$) and \$ ab ) = a a a b aj\§) determinants. The 
subspace of JT^- 1 ) used to solve the PR-EOMCCSD 
problem is spanned by the = | <E>) and |$ J j > ) = 
a b aja,i\§) determinants. By solving Eq. (|35|l . we directly 
obtain the energy differences, w^ 4+1 ' ) = E^ +1 ' — E^ in 
the PA-EOMCCSD case, and w {A ~ l) = E {A ~ l) - E {A) in 
the PR-EOMCCSD case, where E {A+1) and E {A ~ 1] are 
the energies of ground (/i = 0) and excited (fi > 0) states 
of the (A + 1)- and (A — l)-particle systems, respectively, 

and Eq is the ground-state coupled-cluster (in this case, 
CCSD) energy of the A-particle reference system. 

In analogy to the ground-state coupled-cluster theory, 
one can extend the PA-EOMCCSD and PR-EOMCCSD 
schemes to higher excitations. For example, in the nu- 
clear analogs of the quantum-chemical EA-EOMCCSDT 
and IP-EOMCCSDT mHHH methods, which 
we would call here the PA-EOMCCSDT (PA-EOMCC 
singles, doubles, and triples) and PR-EOMCCSDT (PR- 
EOMCC singles, doubles, and triples) approaches, one 
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truncates the ground-state cluster operator at the 
3p-3h clusters T 3 (the M = 3 or CCSDT case) and aug- 
ments the Rff (2p-lh) and R^f 1 \2h-lp) operators, 
Eqs. and respectively, by the 3p-2/i compo- 

nent 



Rn,3p-2h — j2 r Jbc aa(lbaCa k a ji 
in the i?^ 4-1-1 ^ case, and the 3h-2p component 

D 1 ijk c 

-ft/j,3fi-2p = J2 r bc a a a k a j a i 



(36) 



(37) 



in the Rl A 1} case. The resulting PA-EOMCCSDT and 
PR-EOMCCSDT operators R (A+1) and are 

Ri A+1) (3p-2/i) = i^.ip + ^^p-ift + ^ l ,3p-2/ l (38) 

and 

R^ (3h-2 P ) = R^h + R^h-ip + R^3h-2 P , (39) 

respectively. In general, if the cluster operator T^ A > 
defining the ground-state of the A-body system is trun- 
cated at the Mp-Mh component, as shown in Eq. J7|), 
one typically truncates the corresponding PA-EOMCC 
operator and the corresponding PR-EOMCC op- 

erator Rjf~^ as follows: 



M-i 

d(A+1) _ n 
^fl / j I >'iJ,,(n+l)p—nhj 



(40) 



71=0 



where 



R 



1 



fi.(n-\-l)p-nh 



and 



ra!(n+ 1)! 



M-l 
n=0 



up ? 



(41) 



(42) 



where 



R 



fi..(n-\-l)h-np 



n\{n + l)\ 01 -°» 

x a ai ■ ■ ■ a a "ai n ■ ■ ■ a^at, 



(43) 



although, at least in principle, one could extend the 
summations over n in Eqs. (|40|l and (|42|l to X^fiLo 
without affecting the explicit connectedness of the left- 
hand side of the corresponding PA-EOMCC and PR- 
EOMCC eigenvalue problems, Eq. (33). The de- 
tailed discussion of the relationships between trunca- 
tion schemes in the i? p A±1 ' ) and operators in the 
PA-EOMCC or EA-EOMCC and PR-EOMCC or IP- 
EOMCC calculations can be found in Ref. 0] (cf., 



also, Refs. |49j, |90j for additional comments and nu- 
merical tests). The (n + l)p-nh and (n + l)h-np am- 
plitudes faaV.'/.a 1 ,, an d r l a\'.'.' l a n > entering the Rp +1 ^ and 
R^^ operators of Eqs. (f^U|) and l|4*2|> . respectively, 
are obtained by solving the eigenvalue problem defined 
by Eq. i|35|) . in which Hn, open is replaced by Hn{M), 
Eq. (fTOft . The subspace of ( j4+1 ^ relevant to the cor- 
responding truncated PA-EOMCC problem is spanned 



by the |$ a ) and | $°? 1 1 ;" i ° n ) = a a a 



a n°i • • • 



a n a; 



•Oil I*) 



(ri = 1, . . . , Af — 1) determinants. Similarly, the sub- 
space of relevant to the corresponding trun- 
cated PR-EOMCC problem is spanned by the |<I>i) and 

l $ ^::C> = flQ1 ■ • ■ aana ^ ■ ■ ■ «ii«il $ > (n = i, • • • , m-i) 

determinants. 

The PA-EOMCC and PR-EOMCC methods, as de- 
scribed above, and their extensions to two or more va- 
lence particle or holes via multiply attached or multi- 
ply ionized schemes |49j, |76|, |77J, lZS| offer several advan- 
tages compared to the equally accurate, but usually a 
lot more complicated, genuine multi-reference coupled- 
cluster methods of either the valence-universal [9ll l92^| 
or the Hilbert-space or state-universal [9j| type that are 
specifically designed to handle general classes of open- 
shell problems. Although there has been significant 
progress in recent years in the development of gen uine 
multi-reference coup led-cluster theories 0, |9f| |9(| [53, 
IH |H EM HM [lQl DM EM EHE EH > multi-reference 
coupled-cluster calculations are often plagued by intruder 
states, unphysical, singular, and multiple solutions, and 
mathematical difficulties with the proper adaptation of 
the corresponding equations to symmetries of the Hamil- 
tonian if one aims at th e gen eral- purp ose c ompu ter codes 

( C f., e.g., Refs. H OH, Efllliflaliia El ES EUni 

for further information). Some of these issues are cur- 
rentl y being addressed (cf., e.g., Refs. [13, IH, 
llOnl IPOl HM HM HM EHa), but none of these prob- 
lems are present in the PA-EOMCC and PR-EOMCC 
calculations, which could be viewed as the physically 
motivated, intruder-state free, state-selective modifica- 
tions of the powerful and elegant valence-universal multi- 
reference coupled-cluster schemes pioneered by Mukher- 
jee and Lindgren |9lL l92l |. In particular, the use of the 
similarity-transformed Hamiltonian H jv.open obtained in 
the coupled-cluster calculations for an A-particle closed- 
shell system, which commutes with the symmetry oper- 
ators of the original Hamiltonian H, automatically guar- 
antees that the PA-EOMCC and PR-EOMCC wave func- 
tions defined by Eq. (|32|) . obtained by diagonalizing 
Hn \opcn m the appropriate subspaces of Jif^ A+1 > and 
Jf*-" 4 " 1 ), are automatically adapted to the symmetries 
of the Hamiltonian H and one does not have to worry 
about the symmetry-contamination issues facing the im- 
plementations of the open-shell coupled-cluster schemes 
employing the unrestricted reference determinants. 

Our calculations for the ground and low-lying excited 
states of the 15- and 17-particle nuclei around 16 O, re- 
ported in this work, have been performed with the basic 
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PA-EOMCCSD and PR-EOMCCSD methods, in which 
the ground-state of 16 is represented by the CCSD 
wave function e Tl+T2 |$) and the nucleon-attaching and 

nucleon-removing operators Rjf^^ and Rj^ are de- 
fined by Eqs. and $5ty. respectively. We can obtain 
an accurate description of the (A + l)-particle nuclei 17 
and 17 F with the PA-EOMCCSD method, in which we 
include the lp and 2p-lh excitations from the 16 core 
to form the 17-particle systems, since the ground-states 
and the low-lying excited states of the 17 and 17 F nu- 
clei that we have singled out in this work are essentially 
one-quasi-particle states, except, perhaps, for the (3/2)^ 
states of 17 and 17 F, which are resonances. Similarly, 
we can study the (A— l)-particle nuclei 15 and 15 N with 
the basic PR-EOMCCSD approach, in which we include 
the lh and 2h-lp excitations from the 16 closed-shell 
core, since the low-lying states of these nuclei are ex- 
pected to be dominated by one-quasi-hole states with re- 
spect to the A-body reference 16 nucleus. As discussed 
in Ref. 0, there is, for example, almost no experimental 
evidence for the fragmentation of the quasi- hole pi/2 and 
p 3 / 2 states of 16 0. The fact that we use the lp and 2p-\h 
excitations in the PA-EOMCCSD calculations to form 
the (A+ l)-body systems and the fact that we use the lh 
and 2h-lp excitations in the PR-EOMCCSD calculations 
for the (A — l)-body systems mean that we includ e many 
of the same correlations as Fujii et al. |l!5l lll(| . Their 
approach is analo gous to a hermitian coupled-cluster ap- 
proach (see Ref. |116) ). There are, however, differences 
between our PA-EOMCCSD and PR-EOMCCSD calcu- 
lations and the calculations reported by Fujii et al. We 
use a biorthogonal EOMCC formalism, based on diago- 
nalizing the non-hermitian similarity-transformed Hamil- 
tonian f/jv(CCSD), Eq. I|14|) . obtained in highly accurate 
CCSD calculations for the A-body closed-shell nucleus, 
which has been very successful in quantum chemistry and 
molecular physics and which brings a lot of correlations 
within basic truncation schemes, such as EOMCCSD, 
PA-EOMCCSD, and PR-EOMCCSD, through the pres- 
ence of high-order correlation terms in if at (CCSD). We 
also differ in the definition of the model space, since Fu- 
jii et al. use a model space similar to that used in the 
no-core shell-model calculations Q , in which a "triangu- 
lar" energy cutoff is applied to Slater determinants in- 
cluded in the diagonalization of the Hamiltonian, in ad- 
dition to the usual single-particle basis set cutoff. Such 
a model space cannot be used in coupled-cluster calcula- 
tions since it violates the Pauli principle in the summa- 
tions over the intermediate states that emerge through 
products of many-body components of the cluster opera- 
tor in coupled-cluster equations. As mentioned ear- 
lier, the use of a given truncation scheme for the cluster 
operator T^ A > implies specific truncation schemes for the 
EOMCC operators, such as R^ +1) and R^~ X) . Thus, 
we use all lp and 2p-lh or lh and 2h-lp excitations in the 
PA-EOMCCSD and PR-EOMCCSD calculations and all 
Ip-lh and 2p-2h cluster amplitudes t\ and A that are 
allowed by a given single-particle basis set, without im- 



posing additional energy cutoffs on the determinants that 
these excitations correspond to, producing many addi- 
tional and important correlations that are outside model 
spaces used in the no-core shell-model calculations. 

The use of the CCSD method for describing the cor- 
related ground-state of the 16 O reference nucleus in the 
PR-EOMCCSD and PA-EOMCCSD calculations for the 
15- and 17-particle systems around 16 is justified by 
the virtually perfect agreement of the CCSD results and 
the results of the exact shell-model diagonalization of 
the Hamiltonian in the sam e m odel space as used in 
the CCSD calculations [3(1 Jfjjfl- As shown in Refs. 
[H IM HI (cf., also, Ref. [mf), the three-body clus- 
ters T3 play a negligible role in the calculations of the 
binding energy of 16 O when the Hamiltonian includes 
pairwise interactions. They also have almost no effect 
on the lowest-energy J 71 " = 3~ state of 16 O, which is 
dominated by the lp-lh excitations, if the Hamiltonian 
contains the two-body interactions only [371 l39l |40| . In 
this case, the basic EOMCCSD method, defined by Eq. 
in which is replaced by the CCSD operator 

T( A )(CCSD), Eq. Q, and is defined by Eq. JH), 
provides the virtually converged description. Interest- 
ingly enough, our PA-EOMCCSD and PR-EOMCCSD 
results for the binding energies of 17 0/ 17 F and 15 0/ 15 N, 
reported in this work, enable us to explain most of the 
6 MeV difference between the converged coupled-cluster 
and experimental results for the lowest J" = 3~ state of 
16 0, reported, for example, in Ref. [37]], which, according 
to our analysis, may largely be caused by the three-body 
interactions neglected in our calculations that affect the 
relevant single-particle energy spacings (see Sec. lIIIBl for 
a detailed discussion). 

If we were to incorporate T3 components in the ground- 
state coupled-cluster calculations for 16 O, we would also 
have to include the 3p-2h excitations R fJll 3 P -2h, Eq. 13611 . 
in the calculations for the 17-particle nuclei and the 
3h-2p excitations i? /i ,3h-2 P , Eq. 1)3 7|) . in the calcula- 
tions for the 15-particle nuclei. Although we do not 
expect the R^.z^h and i? M ,3/ l -2p terms to be signifi- 
cant for the calculations discussed in this work, we will 
examine their role in the future work, once the corre- 
sponding highly efficient computer codes for large-scale 
nuclear applications are developed. As explained in 
Refs. [i^, [43, it is formally possible to include the 
Rfi,3p-2h and i? Mi 3/ l -2 P components in the PA-EOMCCSD 
and PR-EOMCCSD approaches, which are then called 
the PA-EOMCCSD(3p-2/i) and PR-EOMCCSD(3/i-2p) 
methods, respectively |49j . without losing the explicitly 
connected form of the left-hand side of the corresponding 
PA-EOMCC and PR-EOMCC eigenvalue problems given 
by Eq. i|35|) ■ A promising new development in the area of 
including the relatively expensive R^^p-2h and R^^h-2p 
terms in the PA-EOMCC and PR-EOMCC methods 
is that of the so-called active-space PA-EOMCC (EA- 
EOMCC) and PR-EOMCC (IP-EOMCC) approaches, 
which enable one to reduce the costs of the correspond- 
ing parent PA-EOMCCSDT (or EA-EOMCCSDT), PR- 
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EOMCCSDT (or IP-EOMCCSDT), and other more ex- 
pensive higher-order PA-EOMCC and PR-EOMCC cal- 
culations through a small subset of active single-particle 
states, which in the nuclear physics context would cor- 
respond to the highest occupied and lowest unoccu- 
pied shells of the A-body reference nucleus (see Ref. 
|49| for details). We will investigate the effects of the 
3p-2h and 3h-2p excitations on the PA-EOMCCSD and 
PR-EOMCCSD results for 15 0, 17 0, 15 N, and 17 F re- 
ported in this study using the PA-EOMCCSD (3p-2h), 
PR-EOMCCSD (3h-2p) approaches and the active-space 
variants of the PA-EOMCCSDT and PR-EOMCCSDT 
methods of Ref. |49J in our future works. 

We end this section by presenting the most essential 
algorithmic details of the highly efficient PA-EOMCCSD 
and PR-EOMCCSD computer codes that can be applied 
to valence systems with one valence particle and one va- 
lence hole around the closed shell nucleus, such as 16 O, 
and other computational details pertinent to the specific 
calculations for the 15 0, 17 0, 15 N, and 17 F nuclei dis- 
cussed in this paper. First, as discussed in Sec. Ill Al we 
compute a G matrix in a harmonic oscillator basis. In the 
case of the calculations for the 15-17-particle systems de- 
scribed in this paper, we used model spaces consisting of 
five to eight major oscillator shells; the G matrix and its 
pertinent two-body effective interaction were computed 
for oscillator energies in the range Ml 6 [10, 20] MeV 
and the optimum Ml value for each single-particle ba- 
sis set was determined by finding the minimum on the 
curve representing the dependence of the CCSD ground- 
state energy of 16 O on Ml (see Ref. [35| for the details). 
The renormalized Hamiltonian resulting from the G ma- 
trix calculations was corrected for the spurious center- 
of-mass motion using Eq. @, in which /3coM was cho- 
sen such that the expectation value of the center-of-mass 
Hamiltonian -ffcoM with the ground-state CCSD wave 
function of 16 O is 0.0 MeV. Since the PA-EOMCCSD 
and PR-EOMCCSD methods are based on the linear- 
response-like idea of directly calculating the energy dif- 
ferences ujI A±1) = E {A±1) - E (A) between the (A ± 1)- 
and A-particle systems rather than the total energies 
themselves, where the ground-state of the closed-shell 
A-particle system serves as a reference for the (A ± 1)- 
body systems and where we diagonalize the similarity- 
transformed Hamiltonian obtained in the CCSD calcula- 
tions for the A-particle reference system, we used the 
optimum values of Ml and /3coM determined for 16 O 
in the final PA-EOMCCSD and PR-EOMCCSD calcu- 
lations for 15 0, 17 0, 15 N, and 17 F reported in this pa- 
per. As already mentioned, for larger single-particle ba- 
sis sets with seven or eight major oscillator shells, the 
dependence of the results on Ml is virtually none. The 
dependence of the energies of physical states on 0c,cM is 
virtually none as well (cf. Refs. [M |H S |H EJ for 
the details). Thus, the specific values of Ml and /3coM 
become less and less important when the physical states 
are identified and when larger basis sets are employed (in 
fact, the optimum /3coM values approach zero as the basis 



set increases; cf. Table QJ. Most of our calculations for 
15 0, 17 0, 15 N, and 17 F were based on the N 3 LO nucleon- 
nucleon interaction model of Machleidt and co-workers 
|3(i| , although we also performed the calculations for the 
CD-Bonn interaction model , 28] and the Vis model of the 
Argonne group 27]. As pointed out earlier, the Coulomb 
interaction was included in all of the calculations. 

Once the one- and two-body matrix elements of the 
center-of-mass-corrected effective Hamiltonian, Eq. 
are determined and properly sorted out, as described in 
Sec. Ill Al we set up and solve the CCSD equations for 
the closed-shell A-particle system (in our case, 16 0), the 
PA-EOMCCSD equations for the (A+l)-particle systems 
( 17 and 17 F), and the PR-EOMCCSD equations for the 
(A— l)-particle systems ( 15 and 15 N). Our ground-state 
CCSD computer codes rely on the DIIS solver [74| (see, 
also, Refs. |?3), whereas the PA-EOMCCSD and 
PR-EOMCCSD equations for ground and excited states 
of the (A+ 1)- and (A — l)-particle nucle i are solved with 
the Hirao-Nakatsuji genera lizat ion |l!8| of the Davidson 
diagonalization algorithm 119] to nonhcrmitian eigen- 
value problems. The computationally efficient form of 
the CCSD equations in terms of recursively generated in- 
termediates can be derived diagrammatically using Eqs. 
(|21|l and llL'l'l) . From the point of view of code efficiency, 
it is important to realize that some of the intermediates 
entering the CCSD and other CCSD-based equations rep- 
resent matrix elements of the one- and two-body compo- 
nents of the CCSD similarity-transformed Hamiltonian 
-Hjv,o P on(CCSD), Eq. l|3Tl) . If H n is the n-body compo- 
nent of ^vojfen^' f° r ^ ne one_ an d two-body components 
Hi and H2, respectively, we can write 

H t = hi N[a a a ] (44) 

and 

H 2 = \h^N[a a al 3 a 5 a 1 \, (45) 

where and h?J^ are the one- and two-body matrix ele- 
ments of H^^J that enter the CCSD, PA-EOMCCSD, 
and PR-EOMCCSD equations. As shown in the Ap- 
pendix, matrix elements hP a and h?J^ are calculated using 
the one- and two-body matrix elements of the Hamilto- 
nian in the normal-ordered form, f& and v^fp, respec- 
tively (cf. Eq. [n^i ) and the singly and doubly excited 
cluster amplitudes t l a and i A, defining T± and T2, respec- 
tively. The computationally efficient form of the CCSD 
equations for the case of pairwise interactions in H, in 
terms of selected types of and h 1 ^ and other recur- 
sively generated intermediates is given in the Appendix. 

We also give in the Appendix the computationally effi- 
cient form of the equations defining the PA-EOMCCSD 
and PR-EOMCCSD eigenvalue problems. These are ob- 
tained by applying diagrammatic methods to the PA- 
EOMCCSD and PR-EOMCCSD equations, which can 
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be given the following general form: 

2 



($ a )| [(ffai^,lp)c + X)(flni^,2p-lh)c]|*) 



r ab' 



in the PA-EOMCCSD case, and 



Once the CCSD equations are solved and all one- and 
two-body matrix elements of H^^^ are determined, 
the most expensive steps of the PA-EOMCCSD and PR- 
EOMCCSD methods employing the factorized equations 
shown in the Appendix are n n\ and r? Q v? ul respectively. 

\ > These relatively low, j¥ b — jV & scalings of the costs of 
the CCSD, PA-EOMCCSD, and PR-EOMCCSD calcu- 
lations with the system size (^/K), which are often orders 
of magnitude smaller than the costs of shell-model calcu- 

(47) lations aimed at similar accuracies, are among the most 
important advantages of the coupled-cluster methodol- 
ogy pursued in this work. 



(**\[(H 2 R ll , lh ) c + Y l (B n R li ah-ip)c 

n=l 



]|*> 



(48) 



(49) 



in the PR-EOMCCSD case. Although formally the PA- 
EOMCCSD and PR-EOMCCSD equations require the 
consideration of the three-body components of H^^^ 
(cf. the n = 3 terms in Eqs. lH7|l and JUIJl), we 
do not have to calculate the corresponding six-index 
matrix elements explicitly. With the help of 

diagrammatic techniques, the three-body components 

of H^p^J that enter thc PA-EOMCCSD and PR- 
EOMCCSD equations can be rigorously factorized and 
rewritten in terms of the one- and two-body compo- 

- (CCSD) 

nents of Hj^ n . In consequence, the final working 
equations of' the PA-EOMCCSD and PR-EOMCCSD 
methods in terms of one- and two-body matrix ele- 
ments of the Hamiltonian, f@ and vTg, respectively, 
T\ and T2 cluster amplitudes defining the underlying 
A-particle ground-state CCSD problem, and the -R M ,i p , 
R/j,,2p-ih) R^,ih, and Rfj, t 2h-ip excitation amplitudes defin- 
ing the particle-attaching and particle-removing opera- 



tors, R^ +1 \2p-lh) and R)f~ 1> (2h-lp), respectively, can 
be re-expressed in terms of the one- and two-body matrix 
elements of H^^^ , hP a and h]fp, respectively, and a few 
additional recursively generated intermediates, leading to 
a fully vectorizable algorithm. This computationally ef- 
ficient form of the PA-EOMCCSD and PR-EOMCCSD 
equations is given in the Appendix too. 

In addition to code vectorization, the main advan- 
tage of deriving the CCSD, PA-EOMCCSD, and PR- 
EOMCCSD equations in the form shown in the Appendix 
is the possibility of obtaining the relatively low CPU 
operation count that characterizes these methods. Thc 
CCSD equations and the determination of the full set of 



(A-l), 



one- and two-body matrix elements of H 



(CCSD) 
TV, open 



are char- 



acterized by the n^n^ steps, where, as mentioned ear- 
lier, n Q and n u are the numbers of occupied and unoccu- 
pied orbitals, respectively, in the single-particle basis set. 



III. RESULTS AND DISCUSSION 

A. Results of the PR-EOMCCSD and 
PA-EOMCCSD Calculations with the N 3 LO 
Interaction and their Convergence with the Basis 
Set 

We focus first on the convergence of our PR- 
EOMCCSD and PA-EOMCCSD results for the ground 
and excited states of 15 0, 15 N, 17 and 17 F with the size 
of the single-particle basis set used in the coupled-cluster 
calculations (see Tables HI and ITT|) and address the issue of 
the dependence of these results on the choice of the os- 
cillator parameter Ml. For comparison purposes, we also 
list our previously published ground-state CCSD results 
for 16 O [37j (cf. Table P), since 16 O serves as a refer- 
ence nucleus for thc PR-EOMCCSD and PA-EOMCCSD 
calculations. We limit our discussion of the convergence 
properties of the PR-EOMCCSD and PA-EOMCCSD re- 
sults for 15 0, 15 N, 17 0, and 17 F to the N 3 LO interaction 
model H3. The PR-EOMCCSD and PA-EOMCCSD 
results for the CD-Bonn [2^| and Argonne Vis H3 m ~ 
teractions exhibit almost identical qualitative features in 
terms of their convergence with the number of major os- 
cillator shells and the way they depend on Ml. 

As shown in Tabled the PR-EOMCCSD binding ener- 
gies of 15 and 15 N, the results of the CCSD calculations 
for the binding energy of 16 O, and the PA-EOMCCSD 
binding energies of 17 and 17 F are practically converged 
at the level of eight major oscillator shells. As demon- 
strated earlier for 16 O |3J|, the dependence of these re- 
sults on the oscillator energy Ml 6 [10, 20] MeV is al- 
most negligible, particularly for seven or eight major os- 
cillator shells. In the latter case, the dependence of the 
results on Ml is practically none. This tells us that the 
renormalization of the short-range part of the nucleon- 
nucleon interaction with the no-core G-matrix approach 
combined with the inclusion of singly and doubly excited 
clusters and the corresponding lp, 2p-lh, lh, and 2h-lp 
excitations in thc coupled-cluster and PR-EOMCC/PA- 
EOMCC calculations for thc valence systems around 16 O 
leads to reasonably well converged ground-state energies 
of these systems. It is true that the N 3 LO interaction 
model has a rather soft core, since it carries a cutoff in 
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TABLE I: Total binding energies and binding energies per 
particle (in parentheses) for ls O and 15 N (the PR-EOMCCSD 
values), 16 (the CCSD values), and 17 and 17 F (the PA- 
EOMCCSD values), computed with the N 3 LO interaction 
model 30], as functions of the number of major oscillator 
shells N. All entries (except for the unitless parameter Pcom) 
are m MeV. The results for 16 O are taken from Ref. 1371 . 
The acron ym E xpt stands for the experimental values, taken 
from Ref. |120H . All energies were calculated at the optimum 
values of Ml (the second last row; determined by identifying 
the Ml value at which the CCSD energy of 1 O reaches the 
minimum value) and /3coM (the last row; determined by the 
condition that the expectation value of HcoM with the CCSD 
wave function is 0.0 MeV). For eight major oscillator shells, 
the results are virtually independent of Ml and /3coM = 0.0. 



Nucleus 


N = 5 


N = 6 


N = 7 


N = 8 


Expt 


15 Q 


91.047 


93.772 


93.219 


92.376 


111.955 




(6.070) 


(6.251) 


(6.215) 


(6.158) 


(7.464) 


15 N 


94.127 


96.432 


95.685 


95.086 


115.492 




(6.275) 


(6.429) 


(6.379) 


(6.339) 


(7.699) 


16q 


108.943 


113.341 


112.446 


111.221 


127.619 




(6.810) 


(7.084) 


(7.028) 


(6.951) 


(7.976) 


17 


109.218 


115.534 


115.327 


114.277 


131.762 




(6.425) 


(6.796) 


(6.784) 


(6.722) 


(7.751) 


17 F 


106.097 


112.869 


112.782 


111.510 


128.220 




(6.241) 


(6.640) 


(6.634) 


(6.559) 


(7.542) 


Ml 


13 


11 


10 


11 






1.50 


0.15 


0.05 


0.0 





relative momentum of A = 500 MeV. Thus, in develop- 
ing the effective two-body interaction based on N 3 LO by 
diagonalizing the dcutcron in an oscillator basis, one ob- 
tains a converged result to six leading digits with 50 to 
60 oscillator shells for Ml E [5, 50] MeV. For the CD- 
Bonn and Vig, interactions, one needs more than 100 ma- 
jor shells in order to obtain a converged result for the 
deuteron. However, the advantage of the G-matrix ap- 
proach used in this work is that we can renormalize the 
short-range part of the interaction exactly, since the free 
part of the G matrix is computed in a momentum ba- 
sis first, with the relative momenta |p| G [0, oo). Thus, 
the renormalization problems of the short-range part of 
the two-body interaction, seen, for example, in the no- 
core approach |4|, with a relatively slow convergence as 
a function of the harmonic oscillator excitations, are not 
present here. This means, in turn, that when we use this 
G matrix in coupled-cluster calculations, the results for 
all modern nucleon-nucleon potentials, such as N 3 LO, 
CD-Bonn, and V\% used here, are basically converged 
within eight major shells. It should also be pointed out 
that our re sults for N 3 LO agree very well with those of 
Fujii et al. flT^ . 

The fact that our G-matrix-based coupled-cluster re- 
sults are essentially converged with the basis set has sev- 
eral important consequences for nuclear many-body the- 



ory. We can, for example, claim that any disagreements 
with experimental data are, most likely, due to the miss- 
ing degrees of freedom in our Hamiltonians, such as thrcc- 
nucleon interactions. Indeed, as shown in Table [IJ our 
coupled-cluster calculations miss the experimental bind- 
ing energies (taken from Ref. 120]) by approximately 
1.3 — 1.4 MeV per nucleon for the A = 15 nuclei and by 
approximately 1 MeV per nucleon for the A = 16 and 
A = 17 systems. At least in principle, several factors 
can contribute to these differences, but we believe that 
the three-body interactions are the primary source. It is 
true, for example, that we are using the solution to a two- 
body problem (our G matrix) as the starting point for 
defining a many-body Hamiltonian with pairwise inter- 
actions for the A = 15 — 17 nuclei, and it is known that a 
two-body interaction derived from the diagonalization of 
a three-body problem is different from the corresponding 
two-body interaction derived by diagonali zing the two- 
body problem (e.g., deuteron) P. ll2llll2^ . However, as 
the size of the model space is increased, both two-body 
interactions yield very similar results (see the discussion 
in Ref. |l2lj ). Thus, since we use large model spaces 
with seven or even eight major oscillator shells, the dif- 
ferences between these two types of effective two-body 
interactions are minimal and cannot, as such, contribute 
to the differences between the coupled-cluster and exper- 
imental data observed in our calculations. Clearly, we 
are missing some correlations in our coupled-cluster cal- 
culations, which ignore, for example, T3 clusters in the 
ground-state calculations for 16 and 3p-2h and 3h-2p 
components of and in the PR-EOMCC 

and PA-EOMCC calculations for the 15- and 17-particle 
nuclei, but, as mentioned earlier, the ground and excited 
states of 15 0/ 15 N and 17 0/ 17 F are essentially one-quasi- 
particle and one-quasi-hole states, respectively, in which 
3p-2h and 3h-2 p ex citations play a negligible role. As 
shown in Refs. |37l l3fi IZfT ] . T3 clusters bring in at most 
a total of 1 MeV in the ground-state calculations for 16 O 
(less than 0.1 MeV per nucleon) and cannot, therefore, 
account for the observed differences between the bind- 
ing energies per nucleon. We can thus summarize this 
part of our discussion by stating that the discrepancy 
between experiment and theory observed in Tabled can 
be ascribed to the missing three-body interactions, which 
are not included in our effective Hamiltonians. The ad- 
vantage of the N 3 LO model and similar models based 
on effective field theory, is that they allow for a consis- 
tent deriv ation of three-body terms (see, for example, 
Refs. [3(1 Il23| ). However, we have not developed the 
coupled-cluster codes for dealing with such interactions 
yet. We plan to do it in the future work. 

It is interesting to observe that although the PR- 
EOMCCSD/CCSD/PA-EOMCCSD approaches under- 
bind the five nuclei by about 1-1.4 MeV per particle, 
pointing to the need for the incorporation of three-body 
forces, the relative binding energies of 15 0, 15 N, 16 0, 
17 0, and 17 F obtained in coupled-cluster calculations 
are in good agreement with experiment. For exam- 
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pie, the difference between experimental binding ener- 
gies of 16 and 17 is 0.225 MeV per particle. The 
CCSD and PA-EOMCCSD ground-state energies of 16 
and 17 resulting from the calculations with eight ma- 
jor oscillator shells differ by 0.229 MeV per particle, 
in excellent agreement with experiment. The differ- 
ences between the binding energies for the A = 15 
nuclei and for the A = 17 nuclei are close to the 
experimental values too. From Table [IJ we extract 
BE( 15 N) - BE( 15 0) = 0.181 MeV per particle, when 
the PR-EOMCCSD/N 3 LO approach with eight major 
shells is employed, and 0.235 MeV per particle, when 
the experimental data are used (BE = binding energy). 
Similarly, the PA-EOMCCSD/N 3 LO calculations with 
eight major shells give BE( 17 0) - BE( 17 F) = 0.163 MeV 
per particle, which is in reasonable agreement with the 
experimental result of 0.209 MeV per particle. The 
binding energies per nucleon resulting from the PR- 
EOMCCSD/CCSD/PA-EOMCCSD calculations satisfy 
the ordering 15 < 15 N < 17 F < 17 < 16 0. With 
an exception of the 15 N and 17 F nuclei, whose binding 
energies per particle are close to each other and ordered 
in experiment as 17 F < 15 N, the ordering resulting from 
the PR-EOMCCSD/CCSD/PA-EOMCCSD calculations 
is correct. This clearly illustrates that a great deal of use- 
ful information can be extracted from the relatively in- 
expensive coupled-cluster calculations of the CCSD type 
employing two-body interactions. The total binding en- 
ergies are affected by three-body forces. However, the 
relative binding energies for the nuclei around 16 O seem 
to be reasonably well described at the two-body interac- 
tion level when the coupled-cluster methods are used to 
describe particle correlations. 

We end this subsection by tabulating the results of the 
PR-EOMCCSD and PA-EOMCCSD calculations for the 
low-lying excited states of 15 0, 15 N, 17 0, and 17 F ob- 
tained with the N 3 LO potential (see Table [H] the ex- 
perimental data are taken from Ref. p-24]). Except 
for the (3/2)] 1 " resonance states in 17 and 17 F, all the 
other excited states listed in Table [HI are expected to 
be strongly dominated by one quasi-particle or quasi- 
hole states, meaning that the inclusion of the lp and 
2p-lh excitations in the PA-EOMCCSD calculations and 
the lh and 2h-lp correlations in the PR-EOMCCSD 
calculations should provide a reasonable description of 
these states. This is confirmed in Table [H] The PR- 
EOMCCSD/N 3 LO results for the (3/2)7/ states of 15 O 
and 15 N, employing seven or eight major oscillator shells, 
are particularly impressive, producing errors relative to 
experiment that do not exceed 0.1 MeV. For the parti- 
cle case, the 17 (1/2)+ excited state resulting from the 
PA-EOMCCSD calculations is slightly below the (5/2)+ 
ground state of 17 0, when the N 3 LO interaction and eight 
oscillator shells are employed. Else, the agreement with 
the experimental data is quite satisfactory. Again, as in 
Table U we note a reasonably good convergence in the 
PR-EOMCCSD and PA-EOMCCSD results for the low- 
lying excited states of 15 0, 15 N, 17 0, and 17 F in terms of 



TABLE II: Energies of the low-lying excited states of 15 0, 
15 N, lr O and 17 F, relative to the corresponding ground-state 
energies (the (1/2)7 states of 15 and 15 N and the (5/2)t 
states of lr O and 17 F), computed with the N 3 LO interac- 
tion model 30] and the PR-EOMCCSD ( ls O, 15 N) and PA- 
EOMCCSD ( 17 and 17 F) methods, as functions of the num- 
ber of major oscillator shells TV. All entries are in MeV. Note 
that the experimentally observed (3/2)+ states in lr O and 
17 F are r esonances. The experimental data (Expt) are from 
Ref. |124|| . For the optimum values of hQ, and /3coM , see Table 

m 



Excited state N = 5 N = 6 N = 7 N = 8 Expt 



ls O (3/2)7 


6.515 


6.602 


6.166 


6.264 


6.176 


15 N (3/2)7 


6.354 


6.680 


6.256 


6.318 


6.323 


17 (3/2)+ 


6.298 


6.031 


5.489 


5.675 


5.084 


17 (1/2)+ 


0.328 


0.130 


-0.349 


-0.025 


0.870 


17 F (3/2)+ 


6.460 


6.207 


5.686 


5.891 


5.000 


17 F (1/2)+ 


0.748 


0.544 


0.088 


0.428 


0.495 



the number of harmonic oscillator shells in a basis, which 
is yet another confirmation of the effectiveness of our 
computational procedure, in which we perform coupled- 
cluster calculations with the renormalizcd form of the 
Hamiltonian rather than with the underlying bare inter- 
actions that would lead to a very slow convergence rate 
with the number of single-particle states, making the cal- 
culations unmanageable. 

Within a single-particle picture, the splitting between 
the (3/2)7/ excited and (1/2)7/ ground states in 15 and 
15 N and the splitting between the (3/2)+ excited and 
(5/2)+ ground states in 17 and 17 F should arise from 
the nuclear spin-orbit force. It is interesting to analyze 
to what extent the three-nucleon interactions may affect 
these splittings. The nucleon-nucleon interaction con- 
tains a short-range spin-orbit force, which in a meson- 
exchange model picture originates from heavier vector 
mesons. Several partial waves receive significant contri- 
butions from the two-body spin-orbit force. For example, 
the 3 P2 partial wave, crucial for the pairing properties in 
nuclei and neutron star matter, yields an attractive inter- 
action up to almost 1 GeV in laboratory energy for the 
two-nucleon scattering. This attraction arises from the 
two-body spin-orbit force, since both the central and ten- 
sor force contributions are repulsive. Within the frame- 
work of many-body perturbation theory, the largest con- 
tribution to the spin-orbit force arises from the first-order 
Hartree-Fock diagram. Indeed, for the N3LO model used 
here and for an oscillator energy HVt = 14 MeV, we obtain 
an excitation energy of 5.412 MeV for the 0^3/2 state of 
15 O, in reasonable agreement wit h the experimental and 
coupled-cluster data in Table |H] |125| . At the Hartree- 
Fock diagram level, the origin of the spin-orbit splitting 
comes then from the renormalization of the short range 
two-body spin-orbit force. The nuclear tensor force gives 
also, as a second- and higher-order process, a contribu- 
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tion to the single-particle spin-orbit splitting (see the de- 
tailed discussion in Ref. |2| for further information). The 
authors of Ref. Q show how the second-order diagrams 
in many-body perturbation theory with the 2h-lp and 
2p-lh intermediate states yield repulsive and attractive 
contributions to the single-particle energies, respectively. 
Depending on the strength of the nuclear tensor force, the 
spin-orbit splittings can then be enhanced or reduced. If 
the tensor force is weak, as is the case for the N 3 LO 
model, the reduced higher-order quenching of the tensor 
force terms enhances the spin-orbit splitting with respect 
to the Hartree-Fock diagram. Anticipating the discussion 
in Sec. IIII 01 potentials with a stronger tensor force, such 
as the Vis model of the Argonne group [2]}, result in a 
smaller spin-orbit splitting than the N 3 LO model (and a 
reduction in the (3/2) f - (1/2)7/ and (3/2)+ - (5/2)+ 
spacings in the 15 0/ 15 N and 17 0/ 17 F nuclei, respec- 
tively). The authors of Refs. Q, demonstrated then 
that a two-pion three-nucleon interaction also contributes 
to the spin-orbit splitting. With the inclusion of such 
a term, Pieper and Pandharipande Q], reproduced very 
well the (3/2)^ - (1/2)7/ splitting in 15 N. These find- 
ings were later corroborated by Heisenberg and Mihaila 
in their coupled-cluster calculations with three-body in- 
teractions for 16 (see Refs. (jHl El 0] and the dis " 
cussion in the next subsection as well) . The fact that we 
reproduce very well the experimental (3/2)7/ — (1/2)7/ 
spin-orbit splittings in 15 N and 15 O with the pairwise 
N 3 LO model indicates that the spin-orbit force associ- 
ated with an eventual three-body force for N 3 LO should 
be small. This may be an important finding for our un- 
derstanding of the role of three-body forces in nuclear 
structure calculations. 

In Sec. IIII CI we present results for the binding ener- 
gies and spectra of the 15 N and 15 nuclei and their 17 F 
and 17 O counterparts using the CD-Bonn [2^ and the 
Argonne Vis H3 interaction models as well, so that we 
can see how much the effects due to three-body inter- 
actions may depend on the underlying two-body forces. 
However, before we proceed, let us discuss interesting 
consequences of our PR-EOMCCSD and PA-EOMCCSD 
calculations for 15 N, 15 0, 17 F, and 17 for the nuclear 
structure studies of the excited states of 16 O. 



B. Consequences of the PR-EOMCCSD and 
PA-EOMCCSD Calculations for the Valence 
Systems around 16 O for the Studies of Excitations in 

16q 

Based on the N 3 LO results discussed in the previous 
subsection, we attempt to link our findings to nuclear 
structure studies of the excitations in 16 O. The fact that 
we obtain practically converged results for a given two- 
body Hamiltonian allows us to infer that eventual dis- 
agreements with experiment in the results of ab initio 
calculations for excited states of 16 O can very likely be 
retraced to the degrees of freedom that are not included 



in the existing two-body Hamiltonians. 

Here, we discuss the excited states of 16 with an ex- 
pected lp-lh structure. In Ref. |33|, we obtained con- 
verged results for the lowest-lying 3j~ state of 16 O. For 
the N 3 LO interaction, we obtained an excitation energy 
of about 12 MeV, almost 6 MeV above the experimental 
value of 6.13 MeV. The same Hamiltonian as that used 
here was employed. The low-lying excited states of O 
and, in general, states which involve cross-shell excita- 
tions, have always eluded a proper microscopic descrip- 
tion (see, for example, Refs. p i IHl IB! \im Il3rl Il3l| 
and references therein). 

Let us concentrate on the lowest-energy 3j~ state of 
16 0. In a zero-order approximation, this state may be 
regarded as a state that arises from the single i — > a ex- 
citation from the i = 0pi/ 2 hole state to the a = 0d 5 / 2 
particle state. Relative to the 16 ground state, the en- 
ergy required to produce such an excitation equals 

Ae^ = e v (0d 5/2 ) - e 7r (0p 1/2 ) 

= [BE( 16 0) - BE( 17 F)] + [BE( 16 0) - BE( 15 N)] 
= 11.526 MeV, (50) 

for the proton case, and 

Ae„ = e„(0d 5/2 ) - e v (Qpi/ 2 ) 

= [BE( 16 0) - BE( 17 0)] + [BE( 16 0) - BE( 15 0)] 
= 11.521 MeV, (51) 

for the neutron case, where BE's in the above equations 
represent the relevant total binding energies. In calculat- 
ing the above values of the lp-lh excitation energies Ae w 
and Ae„ that provide us with the zero-order estimates of 
the excitation energy of the lowest 3j~ state of 16 0, we 
used the experimental binding energies listed in Table U 
As we can see from Eqs. 1)500 and <|51|) . the proton and 
neutron excitation energies are practically identical. This 
reflects a well-known feature of the spin-isospin saturated 
systems. Without interactions among nucleons and with 
the above single-particle orbits used as the only active de- 
grees of freedom, all negative parity states with quantum 
numbers J n = 2 _ , 3~ would be at the above energies 
of approximately 11.5 MeV. The interactions among nu- 
cleons lower the energy of the first-excited 3~ state by 
11.5-6.1 = 5.4 MeV. 

Let us now compare the approximate energy spacing 
defining the lowest 3~ state of 16 0, resulting from the 
use of experimental binding energies, as shown above 
(11.5 MeV), with the values of Ae w and Ae„ based on 
the results of coupled-cluster calculations for the bind- 
ing energies of 16 O and valence systems around 16 O ob- 
tained with the N3LO interaction and eight major oscil- 
lator shells. These results are Ae T = 15.846 MeV and 
Ae„ = 15.789 MeV, for proton and neutron excitations, 
respectively, with almost the same difference between the 
proton and neu tron cases as observed in experiment. The 
authors of Ref. [ill obtained 14.72 MeV and 14.64 MeV 
for protons and neutrons, respectively, using the same 
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N 3 LO interaction as used here. We should note, how- 
ever, that their results for 16 are not fully converged 
as a function of the model space size. Using the above 
elementary picture of the Ip-lh excitation defining the 
lowest 3~ state of 16 O, which involves only two orbits 
in the definition of the relevant model space, we can see 
that we are off by approximately 15.8 — 11.5 = 4.3 MeV, 
when we compare the Ae OT and Ae v energy spacings re- 
sulting from coupled-cluster calculations with the exper- 
imental estimates of these spacings. This difference is 
obviously an interaction and method dependent result. 
It is, however, converged as a function of the number 
of oscillator shells in a basis set, showing that the dis- 
crepancy of 4.3 MeV between theory an experiment for 
the energy gap between the Op and IsOd shells accounts 
for a large fraction of the missing 6 MeV needed to re- 
produce the first 3~ state of 16 0. This is, perhaps, the 
most likely candidate for a consistent explanation of the 
large difference between converged coupled-cluster result 
for the lowest 3~ state of 16 and experiment reported 
in Ref. [S^. The above analysis indicates that a large 
fraction of the difference between theory and experiment 
can be traced in this case to errors in reproducing the ex- 
perimental binding energies of 16 O and valence systems 
around 16 O by coupled-cluster methods employing pair- 
wise interactions only. This allows us to conclude that 
a 6 MeV difference between coupled-cluster result and 
experiment for the lowest 3 _ state of 16 is primarily 
caused by the lack of three-body interactions in our cal- 
culations, and much less by the approximate treatment 
of particle correlations by the coupled-cluster methods 
used in our studies. The above analysis also implies that 
with an adjusted gap between the Op and IsOd shells, one 
should be able to get a better reproduction of the excited 
states of 16 which have a well-defined lp-lh structure, 
such as the lowest 3~ state discussed here. One possible 
strategy for describing excited states of closed-shell nu- 
clei dominated by lp-lh excitations might be to keep the 
original two-body Hamiltonian and add additional three- 
body terms via corrections to th e single-pa rticle energies, 
as advocated recently by Zuker |l32l Il33| . 



C. Results for Other Interaction Models 

The binding energies per particle for the three interac- 
tion models examined in this work, namely N 3 LO, CD- 
Bonn, and Vis , are listed in Table IIIII We only show 
the essentially converged results obtained with eight ma- 
jor oscillator shells, since convergence patterns with the 
number of major oscillator shells that characterize the 
N 3 LO, CD-Bonn, and Vis interactions are practically 
identical. 

As expected, the CD-Bonn interaction gives more at- 
traction than N3LO, while the Argonne Vis interaction 
model yields less attraction than the other two models. 
The CD-Bonn potential has the weakest tensor force of 
the three interactions studied here, whereas the Vis in- 



teraction has the strongest tensor force component. It is 
well-known that an interaction model with a weak tensor 
force yields less quenching in the medium for the impor- 
tant 3 Si and 3 Di partial wave contributions to various 
matrix elements of the Hamiltonian. The quenching is as- 
cribed to both a Pauli effect and an energy dependence 
reflected in second- and higher-order terms (see, for ex- 
ample, Ref. |5j for a discussion of this topic in both nuclei 
and nuclear matter). Although all interaction models fit 
properties of the deuteron and the scattering data with a 
X 2 per datum close to 1, the non- localities which are in- 
troduced due to the way the interactions are constructed 
are responsible for different results in a many-body con- 
text. Indeed, the N 3 LO and CD-Bonn models are non- 

TABLE III: A comparison of the binding energies per par- 
ticle for ls O and 15 N (the PR-EOMCCSD values), ie O (the 
CCSD values), and 17 and 17 F (the PA-EOMCCSD values), 
obtained with the N 3 LO I3CS, CD-Bonn |H, and Vis Ha 
potentials, and eight major os cillat or shells, with the exper- 
imental data taken from Ref. 120]. All entries are in MeV. 
For the CD-Bonn and N 3 LO interactions, we used HQ = 11 
MeV. For Vis, we used HQ, = 10 MeV. For eight major shells, 
the results are practically independent of the choice of HQ and 
PcoM = 0.0. 



Interaction 



Nucleus N 3 LO 


CD-Bonn 


Via 


Expt 


15 


6.158 


6.643 


4.789 


7.464 


15 N 


6.339 


6.810 


4.957 


7.699 


16 


6.951 


7.444 


5.469 


7.976 


17 o 


6.722 


7.201 


5.214 


7.751 


17 F 


6.559 


7.048 


5.059 


7.542 



local interactions defined in momentum space. While 
the N 3 LO model is based on chiral Lagrangians with 
nucleons and pions as degrees of freedom, including the 
non-iterative 2ir diagrams at chiral fourth order, the CD- 
Bonn interaction is a traditional meson-exchange model 
that includes the six low-mass mesons n, 5, p,Q,rj and 
the fictitious a meson, which is a 2n resonance. The Ar- 
gonne Vis model is based on a local r-space parametriza- 
tion, dominated by one-pion exchange. The strength of 
the nuclear tensor force is intimately connected with the 
non-localities of the different nucleon-nucleon forces. De- 
pending on how it is quenched in a many-body context, 
one may get less or more attraction. The attractive part 
of, for example, the 3 Si partial wave contribution is more 
attractive in the medium for an interaction with a weak 
tensor than for one with a strong tensor force. Such 
features are clearly seen in the coupled-cluster results re- 
ported in Table ITTll where the potential with the weakest 
tensor force, CD-Bonn, yields more binding than the two 
other models. 

As mentioned previously, our results are practically 
converged as functions of the number of harmonic os- 
cillator shells. Based on our earlier work [33], the triply 
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excited clusters and the related 3p-2h and 3h-2p exci- 
tations in the particle-attaching and particle-removing 

rL A+1) and R { ^ 1] operators of the PA-EOMCC and 
PR-EOMCC theories are expected to have very little im- 
pact on the calculated binding energies. We claim there- 
fore that, except for a small correction due to triples and 
a weak starting energy dependence [3]} , the lack of agree- 
ment between coupled-cluster and experimental binding 
energies is primarily due to the missing physics in our 
Hamiltonians. The main conclusion that one can derive 
from the results of our coupled-cluster calculations with 
different interactions is that every nucleon-nuclcon inter- 
action model needs its own three-body potential. The Ar- 
gonne group has derived sophisticated three-body inter- 
action t erms (see, for example, the extensive elaboration 
of Ref. |134| ~) . The parameters entering their three-body 
interaction models are fitted to reproduce properties of 
light nuclei. These three-body terms follow much of the 
same pion-exchange picture adopted in the construction 
of the Argonne Vis interaction. For the CD-Bonn inter- 
action one would need to derive three-body terms based 
on a meson-exchange p ictur e, as outlined, for example, 
by the Bochum group 135]. However, no such model, 
which accompanies this interaction, has been fully devel- 
oped. The situation for models based on effective field 
theory is much better as three-body terms arise quit e nat - 
urally at given orders in the expansion parameter |l23j| . 
Our coupled-cluster results indicate that every interac- 
tion, due to different non-localities, has its own three- 
body component reflected in different binding energies 
and different spin-orbit splittings (the (3/2)7/ — (1/2)7/ 
spacings in 15 and 15 N and the (3/2)+ - (5/2)+ spac- 
ings in 17 and 17 F), as demonstrated in Table ITTTl which 
lists binding energies per nucleon, and Table IIVI which 
lists the corresponding low-lying excited states of the va- 
lence systems around 16 O examined in this work. 

As shown in Table |IV| the CD-Bonn and the N 3 LO 
models result in the largest spin-orbit splittings (much 
larger than in the case of Vis). I n order to examine this 
behavior in some detail, we have computed all diagrams 
through third order in the G matrix for hil = 14 MeV, 
using many-body perturbation theory as described in 
Ref. |{|, including folded diagrams to infinite order. For 
example, at the Hartree-Fock level, which corresponds to 
the first order in the G matrix, the spin-orbit splittings 
for neutrons between the two hole states in the Op shell 
are 4.85 MeV, 4.41 MeV, and 3.91 MeV for the CD-Bonn, 
N 3 LO and Vis interaction models, respectively. Since 
we are dealing with spin-isospin saturated systems, the 
results for protons are almost the same. The Hartree- 
Fock term yields the largest contribution and receives 
important contributions from the short-range two-body 
spin-orbit force. However, there is also a considerable 
contribution to the splitting that originates from the 
second-order 2h-lp and 2p-lh terms. The correspond- 
ing second-order contributions are 1.81 MeV, 1.73 MeV 
and 1.35 MeV for the same three interactions, respec- 
tively. These perturbation theory estimates agree with 



TABLE IV: A comparison of the energies of the low-lying 
excited states of 15 0, 15 N, 17 and 17 F, relative to the corre- 
sponding ground-state energies (the (1/2) 77 states of ls O and 
15 N and the (5/2)+ states of 17 and 17 F), obtained with 
the PR-EOMCCSD ( ls O and 15 N) and PA-EOMCCSD ( 17 
and 17 F) methods, the N 3 LO CD-Bonn |2||, and Vis E3| 
potentials, and eight major os cillat or shells, with the exper- 
imental data taken from Ref. 124]. All entries are in MeV. 
For the CD-Bonn and N 3 LO interactions, we used Kl = 11 
MeV. For Vis, we used HQ = 10 MeV. For eight major shells, 
the results are practically independent of the choice of HQ and 
/3com = 0.0. 



Interaction 



Excited state N 3 LO 


CD-Bonn 


Vis 


Expt 


15 (3/2)7 


6.264 


7.351 


4.452 


6.176 


15 N (3/2)7/ 


6.318 


7.443 


4.499 


6.323 


17 (3/2)+ 


5.675 


6.406 


3.946 


5.084 


17 (1/2)+ 


-0.025 


0.311 


-0.390 


0.870 


17 F (3/2)+ 


5.891 


6.677 


4.163 


5.000 


17 F (1/2)+ 


0.428 


0.805 


0.062 


0.495 



the ways the (3/2) 1 —(1/2)! spacings in 15 and 15 N and 
the (3/2)+ - (5/2)+ spacings in 17 and 17 F, obtained in 
the corresponding PR-EOMCCSD and PA-EOMCCSD 
calculations, vary with the interaction. This analysis il- 
lustrates, at least to some extent, the role played by the 
quenching of the tensor force via the second and higher- 
order terms in many-body perturbation theory in differ- 
ent interaction models. The perturbative results do not 
stabilize, however, as functions of the oscillator energy, 
a result which is in clos e agr eement with the findings re- 
ported by Fujii et al. |ll5j . With increasing Ml, the 
single-particle splittings increase if one uses an unper- 
turbed harmonic oscillator basis. We defer thus from 
a more elaborate analysis of many-body perturbation 
theory, since it yields results of a rather limited inter- 
est. The problems with many-body perturbation theory, 
such as the lack of a proper indication of convergence in 
terms of G and the difficulties with going beyond third 
order in the interaction, are well known. The coupled- 
cluster methods, including the quantum chemistry in- 
spired CCSD, PA-EOMCCSD, and PR-EOMCCSD ap- 
proximations used in this work, are capable of summing 
large classes of diagrams to infinite order, eliminating 
many of the problems encountered in many-body pertur- 
bation theory calculations, and providing a much more 
stable description of the ground and excited states of 
the valence systems around 16 O. The differences between 
coupled-cluster results obtained with different interac- 
tions point to the need for developing three-body interac- 
tions consistent with a given two-body interaction model. 

It is interesting to note that in spite of the apparent 
differences between the converged coupled-cluster results 
obtained with different pairwise interaction models, the 
relative binding energies of 15 0, 15 N, 16 0, 17 0, and 17 F 
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obtained with different interactions are in good agree- 
ment with experiment and with each other. For example, 
as already mentioned the difference between experimen- 
tal binding energies of 16 and 17 is 0.225 MeV per 
particle. The CCSD and PA-EOMCCSD ground-state 
energies of 16 and 17 resulting from the calculations 
with eight major oscillator shells differ by 0.229 MeV 
per particle for N 3 LO, 0.243 MeV per particle for CD- 
Bonn, and 0.255 MeV per particle for Vis- Similarly, the 
difference between experimental binding energies of 16 O 
and 15 is 0.512 MeV per particle, whereas the CCSD 
and PR-EOMCCSD ground-state energies of 16 and 
15 differ by 0.793, 0.801, and 0.680 MeV per parti- 
cle for the N 3 LO, CD-Bonn, and Vis potentials, respec- 
tively. Here, the differences with experiment are some- 
what greater than in the case of 16 and 17 0, but the 
overall agreement among different potentials is still very 
good. The differences between the binding energies for 
the A = 15 nuclei and for the A = 17 nuclei obtained 
with different interactions are close to one another and 
to the experimental values too. According to Table UTTl 
the experimental value of the binding energy difference 
BE( 15 N) - BE( 15 0) is 0.235 MeV per particle. The PR- 
EOMCCSD calculations with the N 3 LO, CD-Bonn, and 
Vis interactions give 0.181, 0.167, and 0.168 MeV per par- 
ticle, respectively, for the same binding energy difference. 
Similarly, the experimental value of the binding energy 
difference BE( 17 0) - BE( 17 F) is 0.209 MeV per particle. 
The PA-EOMCCSD calculations with the N 3 LO, CD- 
Bonn, and Vis potentials give 0.163, 0.153, and 0.155 
MeV per particle, respectively, for the same binding en- 
ergy difference. In spite of the substantial differences 
between binding energies resulting from the calculations 
with different interactions, which are affected by the 
three-body forces that are expected to be different for 
different pairwise interactions, the binding energies per 
nucleon resulting from our PR-EOMCCSD/CCSD/PA- 
EOMCCSD calculations with eight major oscillator shells 
satisfy 15 O < 15 N < 17 F < 17 < 16 O, independent 
of the interaction used in coupled-cluster calculations. 
This means that once we adjust the value of the bind- 
ing energy of the reference 16 O system, we can obtain 
the interaction-independent ordering of the binding en- 
ergies of the valence nuclei around 16 0. With an ex- 
ception of the 15 N and 17 F nuclei, whose binding en- 
ergy ordering should be reversed, the ordering of binding 
energies per particle resulting from the relatively inex- 
pensive coupled-cluster calculations is in good agreement 
with experiment. These are encouraging findings from 
the point of view of the future applications of coupled- 
cluster methods employing renormalizcd Hamiltonians in 
nuclear physics. 



IV. CONCLUSIONS AND PERSPECTIVES 

We summarize here our main conclusions and perspec- 
tives for future studies. 



1. To our knowledge, this is the first application of 
the ab initio coupled-cluster theory employing the 
renormalized form of the Hamiltonian, combined 
with the PA-EOMCC and PR-EOMCC formalisms 
for open-shell many-fcrmion systems, to nuclear va- 
lence systems with one valence particle or one va- 
lence hole. We have shown that one can obtain 
virtually converged results with given two-body 
Hamiltonians for both binding energies and low- 
lying excited states. This has been possible thanks 
to the development of highly efficient CCSD, PA- 
EOMCCSD, and PR-EOMCCSD computer codes 
and the use of the renormalized Hamiltonians in 
our calculations, which lead to a rapid convergence 
with the number of oscillator shells in a basis. The 
systems whose properties have been studied in this 
work were 15 0, 15 N, 17 and 17 F. An emphasis 
has been placed on states dominated by one-quasi- 
particle configurations. The discrepancies between 
the results of large-scale coupled-cluster calcula- 
tions for these nuclei and the corresponding exper- 
imental data have been traced to the Hamiltonians 
used in the calculations, much less to the correla- 
tions neglected in coupled-cluster approximations 
employed in this study. 

2. Three different nucleon-nucleon interactions have 
been used to define our two-body Hamiltonians. 
These are the N 3 LO model [U, the CD-Bonn in- 
teraction [2^, and the Vis model of the Argonne 
group [23|. All of these interactions yield differ- 
ent binding energies and different energies of the 
excited states. The different binding energies and 
spin-orbit splittings can be related to varying non- 
localities in the nucleon-nucleon interactions. Of 
particular interest here has been the role played by 
the nuclear tensor force. The different behavior of 
the three interaction models examined in this study 
points to the need for the development of the inter- 
action specific three-body forces. 

3. We have also demonstrated that most of the dis- 
crepancy between theory and experiment for the 
lp-lh negative parity states in 16 0, including the 
lowest 3^ state examined in our earlier work |37| . 
can be retraced to the difference between the theo- 
retical and experimental values of the relevant en- 
ergy gaps between neutron or proton states in the 
Op and IsOd shells. 

4. In spite of the differences among interactions, the 
relative binding energies of the 15 0, 15 N, 17 F, 17 0, 
and 16 O resulting from the coupled-cluster calcu- 
lations seem to be virtually independent of the in- 
teraction and in good agreement with experiment. 
The (3/2)7/-(l/2)7/ spacings in 15 and 15 N result- 
ing from the converged coupled-cluster calculations 
with the N 3 LO interaction are in excellent agree- 
ment with experiment, indicating that the spin- 
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orbit force associated with an eventual three-body 
force for N 3 LO should be small. 

There are several obvious extensions to this work. First 
of all, the need for an inclusion of three-body interac- 
tions sets the agenda for forthcoming studies. Moreover, 
it may be useful to examine the role of 3p-2h and 3h-2p 
correlations in the PA-EOMCC and PR-EOMCC calcu- 
lations, which we neglected in this study. For the states 
considered here, the 3p-2h and 3h-2p correlations are ex- 
pected to be small, since the states of 15 0, 15 N, 17 F, and 
17 that we have examined show relatively small depar- 
tures from an independent-particle picture and since the 
underlying T 3 cluster contributions that define the refer- 
ence 16 system are small On the other hand, we 
have tacitly assumed that the 0d 3 / 2 states of 17 F and 17 
are bound states. These states are resonances, and it is 
not yet entirely clear how the non-resonant continuum 
may affect the description of these states. The inclusion 
of such contributions in the description of these states is 
another important point to explore, as demonstrated in 
the recent works on the Gamow sh ell-model and complex- 
scaling techniques |13(1 Il37t Il38| . 
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APPENDIX: FACTORIZED FORM OF THE 
CCSD, PA-EOMCCSD, AND PR-EOMCCSD 
EQUATIONS 

In this appendix, we present the working equations 
defining the CCSD, PA-EOMCCSD, and PR-EOMCCSD 



methods exploited in this study. All of the equations are 
expressed in terms of the one- and two-body matrix el- 
ements of the Hamiltonian in the normal-ordered form, 
fP and uTgi respectively (cf. Eq. (JTTJ; in our case, f£ 

and u^l are the one- and two-body matrix elements of 
the normal-ordered form of the effective Hamiltonian H, 
Eq. J2J), the t\ and t^ b cluster amplitudes defining the 
underlying A-particle ground-state CCSD problem, and, 
in the' case of the PA-EOMCCSD and PR-EOMCCSD 
approaches, the r a , r ab , r l , and r l b amplitudes defin- 
ing the particle-attaching (r a and r a 3 b ) and particle- 
removing (r l and r l b ) operators, i?^ (2p-lh) and 

R ( tl A ~ 1 \2h-lp) , Eqs . (|3ljl and Q^, respectively. As ex- 
plained in Sec imi the CCSD, PA-EOMCCSD, and PR- 
EOMCCSD equations can be cast into a computationally 
efficient factorized form expressed in terms of the one- 
and two-body matrix elements of the CCSD similarity- 
transformed Hamiltonian Hjv )0p en(CCSD), h@ and h^fg, 
respectively, and a few additional intermediates that are 
generated in a recursive manner. The complete set of 
one- and two-body matrix elements of -ff/v i0 pen(CCSD) 
and other intermediates that are needed to set up the 
CCSD, PA-EOMCCSD, and PR-EOMCCSD equations 
is given in Table [V] 

The ground-state CCSD equations for the singly and 
doubly excited cluster amplitudes t l a and t^ b , Eqs. i(2*T|) 
and (EH, can be given the following, computationally 
efficient form: 

Zi — s-i _i_ j'cii Zi tm Ae ±m . Te <mi J_ M e j.'nin 

li a — Ja' 1 a L e n m L a v ma l e > IL m i ea 2 ™i™ ae 

= 0, (A.l) 

T 8 nil e/ T S rati oi A mb b ae 2 m ab J 

= 0. (A.2) 

Note that the left-hand sides of Eqs. lf2*Tl) and (or 
(lA.lll and (|A.2f| N ) represent, respectively, the one- and 
two-body matrix elements h % a and h^ b of Hm (CCSD) (see 
Eqs. (|12fl and (|13I) L The relevant intermediates can be 
found in Table [V] The antisymmetrizers sd vq = £f pq , 
which enter Eq. (|A.2fl and other equations presented in 
this appendix, are defined as 

= ^ = 1 - (pq), (A.3) 

with (pq) representing a transposition of two indices. 
Once the above equations are solved for t\ and t l ^ b , the 
ground-state CCSD energy is calculated using the for- 
mula (cf. Eq. (HU) 

E { A) (M) = {mm + fX + + 2f ati), (A.4) 

which is valid for any truncation scheme M > 2. 

Once the t l a and amplitudes are determined and 
the ground-state CCSD energy of the reference A-hody 
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TABLE V: Explicit algebraic expressions for the one- and 
two-body matrix elements elements of _H"jv i0pen (CCSD) (h a 
and h^fg, respectively) and other intermediates (designated 
by 7) used to construct the computationally efficient form of 
the CCSD, PA-EOMCCSD, and PR-EOMCCSD equations. 

Intermediate Expression™ 



hi 

H 

I be 
'' a t 

h ka 

T cd 
n ab 
Tkl 
' l 'J 

'Ha 

K\ 

Ha 
j,b 

l b 

1 a 

TIC 

1 ab 

ia 

i Jb 

ia 
1 a 

r 



/a 1 Qe in 



je ,m 1 1 ef ,mj 



f3 _|_ v 3 e t m 4 



f"' J 4- h e P 
2 U mi L ef ' ""l l e 



_|_ v ea f k 

cd , 1 cd ,mn 
u ab 1 2 mn b ab 



l~cd im 1 .cd i7n 
^amj'b ~l~ ^bmj'a 



Tjb f m 



rlic j_m 
' ma^b 



' 2 V ij r e/ 
T'jb _ .,cb ,jm 
1 ia u im L ea 

u ab * u ab L e n mb L a T J-rr 

Tc iim 1 Lee iim cc iim ■ 1 Tic ±j\ 

"m^ab ' "bm^ae ^ am^bc 1 2 n 1 " n 

u ia 1 ' L mi b a "to <-e ^ ** IL im b ac 

-i-h e P k + T'3 e f k - !-v ef t jk 

^ i eo ' ia ""e 2 U ai h ef 

f b 4 v be t m 



-v~~ f" 

2 u mn "ea 



v ab * v ab l e i v mb l ae 
fajk _ 1 ujk , 

7) J ° 4- ?; eb 7 li 

2 u mn ' ej 

1 _m.ri 

2 n _ i 



2"mt l » 

l,,eb jm 



"Summation over repeated upper and lower indices is assumed. 
fa = ("1/1/3) and i/^j = (a(3\v\"f8) — {a(3\v\5-y) are the one- 
and two-body matrix elements of the Hamiltonian in the normal- 
ordered form, Eq. 1111 and the t % a and t* J , are the singly and 
doubly excited cluster amplitudes defining the ground-state CCSD 
wave function of the A-body reference system. 



system is known, we can set up and solve the eigen- 
value equations defining the PA-EOMCCSD and PR- 
EOMCCSD methods. The PA-EOMCCSD equations for 



the energy differences J]^ +V) 



= E, 



(A+l) 



-E, 



(A) 



and the 



lp and 2p-lh amplitudes, r a and rj b , respectively, defin- 
ing the ground and excited states of the (A + l)-particle 
system, can be given the following, computationally effi- 
cient, form: 



J 



(<S> a \[H N , opcn (CCSD)Rl A+1 \2p-lh)} c \<S>) - K 



4- h e r m 4. lh ef r 



w p, 1 a 



(A.5) 



($ a }| [ J r7 w , opcn (CCSD)4' 4 + 1 ) (2p-lh)] c \*) = 



^li ' ab- 



LW r m 

1 m ab 



4 n ab' ef 



ma eb 



It f mji 
2 ±mb ab J 



(A.6) 



Similarly, we can use the CCSD values of the singly and 
doubly excited cluster amplitudes defining the ground- 
state wave function of the reference A-body system to 
set up the PR-EOMCCSD eigenvalue equations for the 
energy differences lo^ 1 ' = 



E, 



(A-l) 



Er\ A) and the Ih 



and 2h-lp amplitudes, r % and r lJ b , respectively, defining 
the ground and excited states of the (A — I)-particle 
system. The computationally efficient form of the PR- 
EOMCCSD equations is as follows: 



m [H N , opC n(CCSB)R ( u A -^ (2h-lp)] c \$) 



luie n 
2 ran' 



(A.7) 



21 



(^\[H N , opcn (CCSD)R^\2h-lp)] c m = ^[-\hZ h r™ - V m r m > + \h\r% + |fc«„rT - ^ b r m / + ±7^] 

= 4 A - 1} rt (A.8) 
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